Einstein@Home all-sky search for periodic gravitational waves in LIGO S5 data 



J. Aasi 1 , J. Abadie 1 , B. P. Abbott 1 , R. Abbott 1 , T. D. Abbott 2 , M. Abernathy 3 , T. Accadia 4 , F. Acernese 5a " 5c , 
C. Adams 6 , T. Adams 7 , P. Addesso 8 , R. Adhikari 1 , C. Affeldt 9 ' 10 , M. Agathos lla , K. Agatsuma 12 , P. Ajith 1 , 
B. Allen 9 ' 13 - 10 , A. Allocca 14a ' 14c , E. Amador Ceron 13 , D. Amariutci 15 , S. B. Anderson 1 , W. G. Anderson 13 , 
K. Arai 1 , M. C. Araya 1 , S. Ast 9 - 10 , S. M. Aston 6 , P. Astonc 16a , D. Atkinson 17 , P. Aufmuth 10 ' 9 , C. Aulbcrt 9 ' 10 , 
B. E. Aylott 18 , S. Babak 19 , P. Baker 20 , G. Ballardin 21 , S. Ballmer 22 , Y. Bao 15 , J. C. B. Barayoga 1 , D. Barker 17 , 

F. Baronc 5a ' 5c , B. Barr 3 , L. Barsotti 23 , M. Barsuglia 24 , M. A. Barton 17 , I. Bartos 25 , R. Bassiri 3 ' 26 , M. Bastarrika 3 , 

A. Basti 14a ' 14b , J. Batch 17 , J. Bauchrowitz 9 ' 10 , Th. S. Bauer lla , M. Bebronnc 4 , D. Beck 26 , B. Behnke 19 , 
M. Bejger 27c , M.G. Beker lla , A. S. Bell 3 , C. Bell 3 , I. Belopolski 25 , M. Bcnacquista 28 , J. M. Berliner 17 , 
A. Bcrtolini 9 ' 10 , J. Betzwiescr 6 , N. Beveridge 3 , P. T. Beyersdorf 29 , T. Bhadbade 26 , I. A. Bilenko 30 , G. Billingsley 1 , 
J. Birch 6 , R. Biswas 28 , M. Bitossi 14a , M. A. Bizouard 31a , E. Black 1 , J. K. Blackburn 1 , L. Blackburn 32 , D. Blair 33 , 

B. Bland 17 , M. Blom lla , O. Bock 9 ' 10 , T. P. Bodiya 23 , C. Bogan 9 ' 10 , C. Bond 18 , R. Bondarescu 34 , F. Bondu 35b , 
L. Bonelli 14a < 14b , R. Bonnand 36 , R. Bork 1 , M. Born 9 ' 10 , V. Boschi 14a , S. Bose 37 , L. Bosi 38a , B. Bouhou 24 , 

S. Braccini 14a , C. Bradaschia 14a , P. R. Brady 13 , V. B. Braginsky 30 , M. Branchesi 39a ' 39b , J. E. Brau 40 , J. Breyer 9 ' 10 , 
T. Briant 41 , D. O. Bridges 6 , A. Brillct 35a , M. Brinkmann 9 ' 10 , V. Brisson 31a , M. Britzger 9 ' 10 , A. F. Brooks 1 , 
D. A. Brown 22 , T. Bulik 27b , H. J. Bulten lla ' llb , A. Buonanno 42 , J. Burguet-Castell 43 , D. Buskulic 4 , C. Buy 24 , 
R. L. Byer 26 , L. Cadonati 44 , G. Cagnoli 36 , G. Cagnoli 28 , E. Calloni 5a ' 5b , J. B. Camp 32 , P. Campsie 3 , K. Cannon 45 , 
B. Canuel 21 , J. Cao 46 , C. D. Capano 42 , F. Carbognani 21 , L. Carbone 18 , S. Caride 47 , S. Caudill 48 , M. Cavaglia 49 , 

F. Cavalier 31a , R. Cavalieri 21 , G. Cella 14a , C. Cepeda 1 , E. Cesarini 39b , T. Chalcrmsongsak 1 , P. Charlton 50 , 

E. Chassande-Mottin 24 , W. Chen 46 , X. Chen 33 , Y. Chen 51 , A. Chincarini 52 , A. Chiummo 21 , H. S. Cho 53 , J. Chow 54 

N. Christcnscn 55 , S. S. Y. Chua 54 , C. T. Y. Chung 56 , S. Chung 33 , G. Ciani 15 , F. Clara 17 , D. E. Clark 26 , 
J. A. Clark 44 , J. H. Clayton 13 , F. Cleva 35a , E. Coccia 57a ' 57b , P.-F. Cohadon 41 , C. N. Colacino 14a ' 14b , A. Colla 16a ' 16b 
M. Colombini 16b , A. Conte 16a ' 16b , R. Contc 58 , D. Cook 17 , T. R. Corbitt 23 , M. Cordicr 29 , N. Cornish 20 , A. Corsi 1 , 

C. A. Costa 48 ' 59 , M. Coughlin 55 , J.-P Coulon 35a , P. Couvares 22 , D. M. Coward 33 , M. Cowart 6 , D. C. Coyne 1 , 
J. D. E. Crcighton 13 , T. D. Creighton 28 , A. M. Cruise 18 , A. Cumming 3 , L. Cunningham 3 , E. Cuoco 21 , 

R. M. Cutler 18 , K. Dahl 9 ' 10 , M. Damjanic 9 ' 10 , S. L. Danilishin 33 , S. D'Antonio 57a , K. Danzmann 9 - 10 , V. Dattilo 21 , 
B. Daudcrt 1 , H. Davcloza 28 , M. Davier 31a , E. J. Daw 60 , R. Day 21 , T. Dayanga 37 , R. Dc Rosa 5a ' 5b , D. DeBra 26 , 

G. Dcbreczeni 61 , J. Degallaix 36 , W. Del Pozzo lla , T. Dent 7 , V. Dergachev 1 , R. DeRosa 48 , S. Dhurandhar 62 , 

L. Di Fiore 5a , A. Di Lieto 14a ' 14b , I. Di Palma 9 ' 10 , M. Di Paolo Emilio 57a ' 57c , A. Di Virgilio 14a , M. Diaz 28 , A. Dietz 49 
A. Dietz 4 , F. Donovan 23 , K. L. Dooley 9 ' 10 , S. Doravari 1 , S. Dorsher 63 , M. Drago 64a ' 64b , R. W. P. Drever 65 , 
J. C. Driggers 1 , Z. Du 46 , J.-C. Dumas 33 , S. Dwyer 23 , T. Ebcrlc 9 ' 10 , M. Edgar 3 , M. Edwards 7 , A. Effler 48 , 
P. Ehrens 1 , G. Endroczi 61 , R. Engel 1 , T. Etzel 1 , K. Evans 3 , M. Evans 23 , T. Evans 6 , M. Factourovich 25 , 
V. Fafone 57a - 57b , S. Fairhurst 7 , B. F. Farr 66 , M. Favata 13 , D. Fazi 66 , H. Fchrmann 9 ' 10 , D. Fcldbaum 15 , 
I. Ferrante 14a ' 14b , F. Ferrini 21 , F. Fidecaro 14a ' 14b , L. S. Finn 34 , I. Fiori 21 , R. P. Fisher 22 , R. Flaminio 36 , S. Foley 23 , 
E. Forsi 6 , N. Fotopoulos 1 , J.-D. Fournier 35a , J. Franc 36 , S. Franco 31a , S. Frasca 16a ' 16b , F. Frasconi 14a , M. Frede 9 - 10 , 
M. A. Frei 67 , Z. Frei 68 , A. Freise 18 , R. Frey 40 , T. T. Fricke 9 ' 10 , D. Friedrich 9 ' 10 , P. Fritschel 23 , V. V. Frolov 6 , 
M.-K. Fujimoto 12 , P. J. Fulda 18 , M. Fyffe 6 , J. Gair 69 , M. Galimberti 36 , L. Gammaitoni 38a ' 38b , J. Garcia 17 , 
F. Garufi 5a ' 5b , M. E. Caspar 61 , G. Gelencser 68 , G. Gemme 52 , E. Genin 21 , A. Gennai 14a , L. A. Gergely 70 , 
S. Ghosh 37 , J. A. Giaime 48 ' 6 , S. Giampanis 13 , K. D. Giardina 6 , A. Giazotto 14a , S. Gil-Casanova 43 , C. Gill 3 , 

J. Gleason 15 , E. Goctz 9 ' 10 , G. Gonzalez 48 , M. L. Gorodetsky 30 , S. Gofiler 9 ' 10 , R. Gouaty 4 , C. Graef 9 ' 10 , 
P. B. Graff 69 , M. Granata 36 , A. Grant 3 , C. Gray 17 , R. J. S. Greenhalgh 71 , A. M. Gretarsson 72 , C. Griffo 2 , 
H. Grote 9 ' 10 , K. Grover 18 , S. Gruncwald 19 , G. M. Guidi 39a ' 39b , C. Guido 6 , R. Gupta 62 , E. K. Gustafson 1 , 
R. Gustafson 47 , J. M. Hallam 18 , D. Hammer 13 , G. Hammond 3 , J. Hanks 17 , C. Hanna 1 - 73 , J. Hanson 6 , J. Harms 65 , 
G. M. Harry 74 , I. W. Harry 22 , E. D. Harstad 40 , M. T. Hartman 15 , K. Haughian 3 , K. Hayama 12 , J.-F. Hayau 35b , 
J. Heefner 1 , A. Hcidmann 41 , H. Heitmann 35 , P. Hello 31a , M. A. Hendry 3 , I. S. Heng 3 , A. W. Heptonstall 1 , 
V. Herrera 26 , M. Heurs 9 ' 10 , M. Hewitson 9 - 10 , S. Hild 3 , D. Hoak 44 , K. A. Hodge 1 , K. Holt 6 , M. Holtrop 75 , T. Hong 51 
S. Hooper 33 , J. Hough 3 , E. J. Howell 33 , B. Hughey 13 , S. Husa 43 , S. H. Huttncr 3 , T. Huynh-Dinh 6 , D. R. Ingram 17 , 
R. Inta 54 , T. Isogai 55 , A. Ivanov 1 , K. Izumi 12 , M. Jacobson 1 , E. James 1 , Y. J. Jang 43 , P. Jaranowski 27d , E. Jesse 72 , 
W. W. Johnson 48 , D. I. Jones 76 , R. Jones 3 , R.J.G. Jonker lla , L. Ju 33 , P. Kalmus 1 , V. Kalogera 66 , S. Kandhasamy 63 

G. Kang 77 , J. B. Kanner 42 ' 32 , M. Kasprzack 21 ' 31a , R. Kasturi 78 , E. Katsavounidis 23 , W. Katzman 6 , H. Kaufer 9 ' 10 , 
K. Kaufman 51 , K. Kawabe 17 , S. Kawamura 12 , F. Kawazoe 9 ' 10 , D. Keitel 9 ' 10 , D. Kelley 22 , W. Kells 1 , D. G. Keppel 1 

Z. Keresztes 70 , A. Khalaidovski 9 ' 10 , F. Y. Khalili 30 , E. A. Khazanov 79 , B. K. Kim 77 , C. Kim 80 , H. Kim 9 ' 10 , 
K. Kim 81 , N. Kim 26 , Y. M. Kim 53 , P. J. King 1 , D. L. Kinzel 6 , J. S. Kissel 23 , S. Klimenko 15 , J. Kline 13 , 



2 



K. Kokeyama 48 , V. Kondrashov 1 , S. Koranda 13 , W. Z. Korth 1 , I. Kowalska 27b , D. Kozak 1 , V. Kringel 9 - 10 , 
B. Krishnan 19 , A. Kr61ak 27a > 27c , G. Kuchn 9 ' 10 , P. Kumar 22 , R. Kumar 3 , R. Kurdyumov 26 , P. Kwee 23 , P. K. Lam 54 , 
M. Landry 17 , A. Langley 65 , B. Lantz 26 , N. Lastzka 9 * 10 , C. Lawrie 3 , A. Lazzarini 1 , P. Lead 19 , C. H. Lee 53 , 



46 



H. K. Lee 81 , H. M. Lee 82 , J. R. Leong 9 ' 10 , I. Leonor 40 , N. Leroy 31a , N. Letendre 4 , V. Lhuillier 17 , J. Li 
T. G. F. Li lla , P. E. Lindquist 1 , V. Litvine 1 , Y. Liu 46 , Z. Liu 15 , N. A. Lockerbie 83 , D. Lodhia 18 , J. Logue 3 , 
M. Lorcnzini 39a , V. Lorictte 31b , M. Lormand 6 , G. Losurdo 39a , J. Lough 22 , M. Lubinski 17 , H. Luck 9 * 10 , 

A. P. Lundgren 9 ' 10 , J. Macarthur 3 , E. Macdonald 3 , B. Machenschalk 9 ' 10 , M. Maclnnis 23 , D. M. Macleod 7 , 
M. Mageswaran 1 , K. Mailand 1 , E. Majorana 16a , I. Maksimovic 31b , V. Malvezzi 57a , N. Man 35a , I. Mandel 18 , 
V. Mandic 63 , M. Mantovani 14a < 14c , F. Marchesoni 38a , F. Marion 4 , S. Marka 25 , Z. Marka 25 , A. Markosyan 26 , 

E. Maros 1 , J. Marque 21 , F. Martelli 39a ' 39b , I. W. Martin 3 , R. M. Martin 15 , J. N. Marx 1 , K. Mason 23 , A. Masserot 4 , 

F. Matichard 23 , L. Matone 25 , R. A. Matzncr 84 , N. Mavalvala 23 , G. Mazzolo 9 ' 10 , R. McCarthy 17 , D. E. McClelland 54 , 

S. C. McGuire 85 , G. Mclntyre 1 , J. Mclver 44 , G. D. Meadors 47 , M. Mehmet 9 ' 10 , T. Meier 10 ' 9 , A. Mclatos 56 , 
A. C. Melissinos 86 , G. Mendell 17 , D. F. Menendez 34 , R. A. Mercer 13 , S. Meshkov 1 , C. Messenger 7 , M. S. Meyer 6 , 
H. Miao 51 , C. Michel 36 , L. Milano 5a ' 5b , J. Miller 54 , Y. Minenkov 57a , C. M. F. Mingarclli 18 , V. P. Mitrofanov 30 , 
G. Mitselmakher 15 , R. Mittleman 23 , B. Moe 13 , M. Mohan 21 , S. R. P. Mohapatra 44 , D. Moraru 17 , G. Moreno 17 , 

N. Morgado 36 , A. Morgia 57a < 57b , T. Mori 12 , S. R. Morriss 28 , S. Mosca 5a ' 5b , K. Mossavi 9 ' 10 , B. Mours 4 , 
C. M. Mow-Lowry 54 , C. L. Mueller 15 , G. Mueller 15 , S. Mukherjee 28 , A. Mullavey 48 ' 54 , H. Miiller-Ebhardt 9 ' 10 , 
J. Munch 87 , D. Murphy 25 , P. G. Murray 3 , A. Mytidis 15 , T. Nash 1 , L. Naticchioni 16a < 16b , V. Necula 15 , J. Nelson 3 , 
I. Neri 38a ' 38b , G. Newton 3 , T. Nguyen 54 , A. Nishizawa 12 , A. Nitz 22 , F. Nocera 21 , D. Nolting 6 , M. E. Normandin 28 , 
L. Nuttall 7 , E. Ochsner 13 , J. O'Dell 71 , E. Oelker 23 , G. H. Ogin 1 , J. J. Oh 88 , S. H. Oh 88 , R. G. Oldenberg 13 , 
B. O'Reilly 6 , R. O'Shaughnessy 13 , C. Osthelder 1 , C. D. Ott 51 , D. J. Ottaway 87 , R. S. Ottens 15 , H. Ovcrmier 6 , 

B. J. Owen 34 , A. Page 18 , L. Palladino 57a > 57c , C. Palomba 16a , Y. Pan 42 , F. Paoletti 14a ' 21 , R. Paoletti 14a , 
M. A. Papa 19 ' 13 , M. Parisi 5a < 5b , A. Pasqualetti 21 , R. Passaquieti 14a < 14b , D. Passuello 14a , M. Pedraza 1 , S. Pcnn 78 , 
A. Perreca 22 , G. Persichctti 5a < 5b , M. Phelps 1 , M. Pichot 35a , M. Pickenpack 9 - 10 , F. Piergiovanni 39a < 39b , V. Pierro 8 , 
M. Pihlaja 63 , L. Pinard 36 , I. M. Pinto 8 , M. Pitkin 3 , H. J. Pletsch 9 * 10 , M. V. Plissi 3 , R. Poggiani 14a < 14b , J. P61d 9 - 10 , 
F. Postiglione 58 , C. Poux 1 , M. Prato 52 , V. Prcdoi 7 , T. Prcstcgard 63 , L. R. Price 1 , M. Prijatelj 9 - 10 , M. Principe 8 , 
S. Privitera 1 , R. Prix 9 - 10 , G. A. Prodi 64a - 64b , L. G. Prokhorov 30 , O. Puncken 9 ' 10 , M. Punturo 38a , P. Puppo 16a , 
V. Quctschkc 28 , R. Quitzow- James 40 , F. J. Raab 17 , D. S. Rabeling lla < llb , I. Racz 61 , H. Radkins 17 , P. Raffai 25 * 68 , 
M. Rakhmanov 28 , C. Ramct 6 , B. Rankins 49 , P. Rapagnani 16a > 16b , V. Raymond 66 , V. Re 57a ' 57b , C. M. Reed 17 , 

T. Reed 89 , T. Regimbau 35a , S. Rcid 3 , D. H. Reitze 1 , F. Ricci 16a > 16b , R. Riesen 6 , K. Riles 47 , M. Roberts 26 , 
N. A. Robertson 1 - 3 , F. Robinet 31a , C. Robinson 7 , E. L. Robinson 19 , A. Rocchi 57a , S. Roddy 6 , C. Rodriguez 66 , 
M. Rodruck 17 , L. Rolland 4 , J. G. Rollins 1 , J. D. Romano 28 , R. Romano 5a - 5c , J. H. Romic 6 , D. Rosihska 27c < 27f , 
C. Rover 9 ' 10 , S. Rowan 3 , A. Rudiger 9 ' 10 , P. Ruggi 21 , K. Ryan 17 , F. Salemi 9 - 10 , L. Sammut 56 , V. Sandberg 17 , 
S. Sankar 23 , V. Sannibale 1 , L. Santamarfa 1 , I. Santiago-Prieto 3 , G. Santostasi 90 , E. Saracco 36 , B. S. Sathyaprakash 7 , 
P. R. Saulson 22 , R. L. Savage 17 , R. Schilling 9 ' 10 , R. Schnabel 9 ' 10 , R. M. S. Schofield 40 , B. Schulz 9 - 10 , 
B. F. Schutz 19 ' 7 , P. Schwinberg 17 , J. Scott 3 , S. M. Scott 54 , F. Seifert 1 , D. Sellers 6 , D. Sentenac 21 , A. Sergeev 79 , 
D. A. Shaddock 54 , M. Shaltev 9 ' 10 , B. Shapiro 23 , P. Shawhan 42 , D. H. Shoemaker 23 , T. L Sidery 18 , X. Siemens 13 , 
D. Sigg 17 , D. Simakov 9 - 10 , A. Singer 1 , L. Singer 1 , A. M. Sintcs 43 , G. R. Skelton 13 , B. J. J. Slagmolcn 54 , 
J. Slutsky 48 , J. R. Smith 2 , M. R. Smith 1 , R. J. E. Smith 18 , N. D. Smith-Lefebvre 15 , K. Somiya 51 , B. Sorazu 3 , 
F. C. Speirits 3 , L. Sperandio 57a < 57b , M. Stefszky 54 , E. Steinert 17 , J. Steinlechner 9 - 10 , S. Steinlechner 9 ' 10 , 
S. Steplewski 37 , A. Stochino 1 , R. Stone 28 , K. A. Strain 3 , S. E. Strigin 30 , A. S. Stroeer 28 , R. Sturani 39a < 39b , 
A. L. Stuvcr 6 , T. Z. Summcrscales 91 , M. Sung 48 , S. Susmithan 33 , P. J. Sutton 7 , B. Swinkels 21 , G. Szeifert 68 , 
M. Tacca 21 , L. Taffarello 64c , D. Talukder 37 , D. B. Tanner 15 , S. P. Tarabrin 9 ' 10 , R. Taylor 1 , A. P. M. ter Braack lla , 
P. Thomas 17 , K. A. Thorne 6 , K. S. Thorne 51 , E. Thrane 63 , A. Thuring 10 ' 9 , C. Titsler 34 , K. V. Tokmakov 83 , 

C. Tomlinson 60 , A. Toncelli 14a > 14b , M. Tonclli 14a < 14b , O. Torre 14a < 14c , C. V. Torres 28 , C. I. Torrie 1 ' 3 , 
E. Tourncficr 4 , F. Travasso 38a < 38b , G. Traylor 6 , M. Tse 25 , D. Ugolini 92 , H. Vahlbruch 10 - 9 , G. Vajente 14a < 14b , 
J. F. J. van den Brand lla ' llb , C. Van Den Broeck lla , S. van der Putten lla , A. A. van Veggel 3 , S. Vass 1 , 
M. Vasuth 61 , R. Vaulin 23 , M. Vavoulidis 31a , A. Vecchio 18 , G. Vcdovato 64c , J. Vcitch 7 , P. J. Veitch 87 , 
K. Venkateswara 93 , D. Verkindt 4 , F. Vctrano 39a > 39b , A. Vicere 39a < 39b , A. E. Villar 1 , J.-Y. Vinct 35a , S. Vitale lla , 
H. Vocca 38a , C. Vorvick 17 , S. P. Vyatchanin 30 , A. Wade 54 , L. Wade 13 , M. Wade 13 , S. J. Waldman 23 , L. Wallace 1 , 
Y. Wan 46 , M. Wang 18 , X. Wang 46 , A. Wanner 9 - 10 , R. L. Ward 24 , M. Was 31a , M. Weinert 9 ' 10 , A. J. Weinstein 1 , 
R. Weiss 23 , T. Welborn 6 , L. Wen 51 < 33 , P. Wessels 9 ' 10 , M. West 22 , T. Westphal 9 ' 10 , K. Wette 9 - 10 , J. T. Whelan 67 , 
S. E. Whitcomb 1 - 33 , D. J. White 60 , B. F. Whiting 15 , K. Wiesner 9 - 10 , C. Wilkinson 17 , P. A. Willems 1 , 



3 



L. Williams 15 , R. Williams 1 , B. Willke 9 ' 10 , M. Wimmer 9 ' 10 , L. Winkelmann 9 ' 10 , W. Winkler 9 ' 10 , C. C. Wipf 23 , 
A. G. Wiseman 13 , H. Witter 9 ' 10 , G. Woan 3 , R. Wooley 6 , J. Worden 17 , J. Yablon 66 , I. Yakushin 6 , H. Yamamoto 1 , 
K. Yamamoto 64b ' 64d , C. C. Yancey 42 , H. Yang 51 , D. Yeaton-Massey 1 , S. Yoshida 94 , M. Yvert 4 , A. Zadrozny 27c , 
M. Zanolin 72 , J.-P. Zendri 64c , F. Zhang 46 , L. Zhang 1 , C. Zhao 33 , N. Zotov 89 , M. E. Zuckcr 23 , and J. Zwcizig 1 
(The LIGO Scientific Collaboration and the Virgo Collaboration) 

and D. P. Anderson 95 

1 LIGO - California Institute of Technology, Pasadena, CA 91125, USA 
2 California State University Fullerton, Fullerton CA 92831 USA 
3 SUPA, University of Glasgow, Glasgow, G12 8QQ, United Kingdom 
4 Laboratoire d'Annecy-le-Vieux de Physique des Particules (LAPP), 
Universite de Savoie, CNRS/IN2P3, F-74941 Annecy-Le-Vieux, France 
5 INFN, Sezione di Napoli a ; Universita di Napoli Federico 11^ Complesso Universitario di Monte S.Angelo, 
1-80126 Napoli; Universita di Salerno, Fisciano, 1-84084 Salerno c , Italy 
6 LIGO - Livingston Observatory, Livingston, LA 70754, USA 
7 Cardiff University, Cardiff, CF24 3 A A, United Kingdom 
8 University of Sannio at Benevento, 1-82100 Benevento, Italy and INFN (Sezione di Napoli), Italy 
9 Albert- Einstein- Institut, Max- Planck- Institut fur Gravitationsphysik, D-30167 Hannover, Germany 
10 Leibniz Universitdt Hannover, D-30167 Hannover, Germany 
11 Nikhef Science Park, Amsterdam, the Netherlands 0, ; VU University Amsterdam, 
De Boelelaan 1081, 1081 HV Amsterdam, the Netherlands b 
12 National Astronomical Observatory of Japan, Tokyo 181-8588, Japan 
13 University of Wisconsin-Milwaukee, Milwaukee, WI 53201, USA 
14 INFN, Sezione di Ptsa a ; Universita di Pisa b ; 1-56127 Pisa; Universita di Siena, 1-53100 Siena c , Italy 
15 University of Florida, Gainesville, FL 32611, USA 
16 INFN, Sezione di Roma a ; Universita 'La Sapienza* 1 , 1-00185 Roma, Italy 
17 LIGO - Hanford Observatory, Richland, WA 99352, USA 
18 University of Birmingham, Birmingham, B15 2TT, United Kingdom 
19 Albert- Einstein- Institut, Max- Planck- Institut fur Gravitationsphysik, D-14476 Golm, Germany 
20 Montana State University, Bozeman, MT 59717, USA 
21 European Gravitational Observatory (EGO), 1-56021 Cascina (PI), Italy 
22 Syracuse University, Syracuse, NY 13244, USA 
23 LIGO - Massachusetts Institute of Technology, Cambridge, MA 02139, USA 
24 APC, AstroParticule et Cosmologie, Universite Paris Diderot, 
CNRS/IN2P3, CEA/Irfu, Observatoire de Paris, Sorbonne Paris Cite, 
10, rue Alice Domon et Leonie Duquet, 75205 Paris Cedex 13, France 
25 Columbia University, New York, NY 10027, USA 
26 Stanford University, Stanford, CA 94305, USA 
27 IM-PAN 00-956 Warsaw a ; Astronomical Observatory^ Warsaw University 00-478 Warsaw ; CAMK-PAN 00-716 Warsaw ; 
Bialystok University 15-424 Bialystok d ; NCBJ 05-400 Swierk-Otwock e ; Institute of Astronomy 65-265 Zielona Gora' , Poland 
28 The University of Texas at Brownsville, Brownsville, TX 78520, USA 
29 San Jose State University, San Jose, CA 95192, USA 
30 'Moscow State University, Moscow, 119992, Russia 
31 LAL, Universite Paris-Sud, IN2P3/CNRS, F-91898 Orsay° ; ESPCI, CNRS, F-75005 Paris , France 
32 NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA 
University of Western Australia, Crawley, WA 6009, Australia 
34 The Pennsylvania State University, University Park, PA 16802, USA 
3S 'Universite Nice-Sophia- Antipolis, CNRS, Observatoire de la Cote d'Azur, 
F-06304 Nice ; Institut de Physique de Rennes, CNRS, Universite de Rennes 1, 35042 Rennes b , France 
36 Laboratoire des Materiaux Avances (LMA), IN2P3/CNRS, F-69622 Villeurbanne, Lyon, France 
37 Washington State University, Pullman, WA 99164, USA 
3S INFN, Sezione di Perugia" ; Universita di Perugia , 1-06123 Perugia, Italy 
39 INFN, Sezione di Firenze, 1-50019 Sesto Fiorentino a ; Universita degli Studi di Urbino 'Carlo Bo', 1-61029 Urbmo b , Italy 

40 University of Oregon, Eugene, OR 97403, USA 
41 Laboratoire Kastler Brossel, ENS, CNRS, UPMC, 
Universite Pierre et Marie Curie, 4 Place Jussieu, F-75005 Paris, France 
42 University of Maryland, College Park, MD 20742 USA 

43 Universitat de les Illes Balears, E-07122 Palma de Mallorca, Spain 

44 University of Massachusetts - Amherst, Amherst, MA 01003, USA 

45 Canadian Institute for Theoretical Astrophysics, 
University of Toronto, Toronto, Ontario, M5S 3H8, Canada 
46 Tsinghua University, Beijing 100084 China 



4 



47 University of Michigan, Ann Arbor, MI 48109, USA 
48 Louisiana State University, Baton Rouge, LA 70803, USA 
49 The University of Mississippi, University, MS 38677, USA 
50 Charles Sturt University, Wagga Wagga, NSW 2678, Australia 
51 Caltech-CaRT, Pasadena, CA 91125, USA 
52 INFN, Sezione di Genova; 1-16146 Genova, Raly 
53 Pusan National University, Busan 609-735, Korea 
54 Australian National University, Canberra, ACT 0200, Australia 
55 Carleton College, Northfield, MN 55057, USA 
56 The University of Melbourne, Parkville, VIC 3010, Australia 
57 INFN, Sezione di Roma Tor Vergata a ; Universita di Roma Tor Vergata, 
1-00133 Roma b ; Universita dell'Aquila, 1-67100 L'Aquila c , Italy 
58 University of Salerno, 1-84084 Fisciano (Salerno), Italy and INFN (Sezione di Napoli), Italy 
59 Instituto Nacional de Pesquisas Espaciais, 12227-010 - Sao Jose dos Campos, SP, Brazil 
60 The University of Sheffield, Sheffield S10 2TN, United Kingdom 
61 WIGNER RCP, RMKI, H-1121 Budapest, Konkoly Thege Miklos ut 29-33, Hungary 
62 Inter- University Centre for Astronomy and Astrophysics, Pune - 4H007, India 
63 University of Minnesota, Minneapolis, MN 55455, USA 
64 INFN, Gruppo Collegato di Trento" and Universita di Trento b , 1-38050 Povo, Trento, 
Italy; INFN, Sezione di Padova c and Universita di Padova d , 1-35131 Padova, Italy 
65 California Institute of Technology, Pasadena, CA 91125, USA 

66 Northwestern University, Evanston, IL 60208, USA 
67 Rochester Institute of Technology, Rochester, NY 14623, USA 
68 Edtvds Lor and University, Budapest, 1117 Hungary 
69 University of Cambridge, Cambridge, CB2 1 TN, United Kingdom 
70 University of Szeged, 6720 Szeged, Dom ter 9, Hungary 
71 Rutherford Appleton Laboratory, HSIC, Chilton, Didcot, Oxon 0X11 OQX United Kingdom 
72 Embry- Riddle Aeronautical University, Prescott, AZ 86301 USA 
73 Perimeter Institute for Theoretical Physics, Ontario, N2L 2Y5, Canada 
74 American University, Washington, DC 20016, USA 
75 University of New Hampshire, Durham, NH 03824, USA 
76 University of Southampton, Southampton, S017 IB J, United Kingdom 
77 Korea Institute of Science and Technology Information, Daejeon 305-806, Korea 
78 Hobart and William Smith Colleges, Geneva, NY 14456, USA 
79 Institute of Applied Physics, Nizhny Novgorod, 603950, Russia 
80 Lund Observatory, Box 43, SE-221 00, Lund, Sweden 
81 Hanyang University, Seoul 133-791, Korea 
82 Seoul National University, Seoul 151-742, Korea 

83 University of Strathclyde, Glasgow, Gl 1XQ, United Kingdom 

84 The University of Texas at Austin, Austin, TX 78712, USA 

85 Southern University and A&M College, Baton Rouge, LA 70813, USA 

86 University of Rochester, Rochester, NY 14627, USA 

87 University of Adelaide, Adelaide, SA 5005, Australia 

88 National Institute for Mathematical Sciences, Daejeon 305-390, Korea 
89 Louisiana Tech University, Ruston, LA 71272, USA 
90 McNeese State University, Lake Charles, LA 70609 USA 
91 Andrews University, Berrien Springs, MI 49104 USA 
92 Trinity University, San Antonio, TX 78212, USA 
93 University of Washington, Seattle, WA, 98195-4290, USA 
94 Southeastern Louisiana University, Hammond, LA 70402, USA 
95 University of California at Berkeley, Berkeley, CA 94720 USA 
(Dated: August 7, 2012) 

This paper presents results of an all-sky search for periodic gravitational waves in the frequency 
range [50, 1 190] Hz and with frequency derivative range of ~ [-20, 1.1] x 10~ 10 Hz s" 1 for the 
fifth LIGO science run (S5). The search uses a non-coherent Hough-transform method to combine 
the information from coherent searches on timescales of about one day. Because these searches 
are very computationally intensive, they have been carried out with the Einstein@Home volunteer 
distributed computing project. Post-processing identifies eight candidate signals; deeper follow-up 
studies rule them out. Hence, since no gravitational wave signals have been found, we report upper 
limits on the intrinsic gravitational wave strain amplitude ho- For example, in the 0.5 Hz-wide band 
at 152.5 Hz, we can exclude the presence of signals with ho greater than 7.6 x 10~ 25 at a 90% 
confidence level. This search is about a factor 3 more sensitive than the previous EinsteinQHome 
search of early S5 LIGO data. 
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I. INTRODUCTION 

A promising class of sources for detectable gravita- 
tional wave signals are rapidly rotating neutron stars 
with non-axisymmetric deformations [TH5]. Such objects 
are expected to emit long-lived continuous-wave (CW) 
signals. In the rest frame of the neutron star, these waves 
have a constant amplitude and are quasi-monochromatic 
with a slowly decreasing intrinsic frequency. They are 
received at Earth-based detectors with a Doppler modu- 
lation due to the relative motion between the source and 
the detector. Consequently the observed phase evolution 
depends on the intrinsic signal frequency, the first fre- 
quency time-derivative (also called spindown), and the 
source sky position; these parameters shall collectively 
be called the phase evolution parameters. While using 
higher order frequency derivatives could potentially im- 
prove the astrophysical detection efficiency in a part of 
the parameter space (see Sec. Ill I, wc shall not consider 
them in this paper for computational reasons. Finally, 
the received signal has a time-dependent amplitude mod- 
ulation due to the (time-dependent) relative geometry of 
the wave and the detector. 

The previous two decades have seen the construction 
and operation of several kilometer-scale laser interfero- 
metric gravitational wave detectors [6HH] . The detectors 
and the data analysis tools have steadily improved over 
this period. These have made it possible to search for 
various gravitational wave signals with ever-improving 
sensitivity. In this paper we focus on data from the fifth 
science run (S5) of the LIGO (Laser Interferometer Grav- 
itational wave Observatory) detectors, collected between 
the GPS times of 815 155 213 s (Fri Nov 04 16:00:00 UTC 
2005) and 875 145 614 s (Sun Sep 30 00:00:00 UTC 2007). 
The LIGO network [B] consists of three detectors: one at 
Livingston, Louisiana, USA, with an arm length of 4 km 
(L) and two in the same vacuum envelope at Hanford, 
Washington, USA, with arm lengths of 4 km (H) and 
2 km, respectively. Only data from H and L detectors 
are used in this paper. The Virgo and GEO 600 detec- 
tors also collected data during the same time interval but 
are not used in this analysis, which is optimized for two 
detectors with similar sensitivities. 

A coherent strategy for extracting faint CW signals 
buried in noisy data using standard maximum-likelihood 
techniques in the presence of "nuisance parameters" was 
derived in [12"] . The resulting detection statistic is the 
so-called ^-statistic, which has since been generalized to 
the case of multiple detectors [HI [14] . The ^-statistic 
has also been shown to arise as a special case in a more 
general Bayesian framework [15] , Using the ^-statistic 
means that we need to search explicitly only over the 
phase evolution parameters. 

Coherent wide-parameter-space searches utilizing the 
^-statistic have been carried out since the second LIGO 
science run JT6J El]- The amplitude sensitivity of this 
type of search improves as the square root of the time 
baseline. However, the template bank spacing decreases 



dramatically with the time baseline, and the computa- 
tional requirements of the search increase rapidly. Even 
with a coherent time baseline of just a few days, a 
widc-frequency-band all-sky search is computationally 
extremely challenging. It becomes completely unfeasi- 
ble if one considers instead time baselines on the order 
of months. 

As is often the case with computationally bound prob- 
lems, hierarchical approaches have been proposed [TB1 - 
|2"0] . In these strategies, the entire data set is split into 
shorter segments. Each segment is analyzed coherently, 
and afterwards the information from the different seg- 
ments is combined incoherently (which means that the 
phase information is lost). The amplitude sensitivity 
grows at best with the fourth root of the number of 
segments. Such methods have been used in previous 
wide-parameter-space searches published by the LIGO 
and Virgo collaborations [2TH26] . 

A subset of these searches [23 HH [551 121] used seg- 
ments sufficiently short (1800 s) that the signal remains 
within a single Fourier frequency bin in each segment. 
In these simple Fourier transform suffices for 

each segment. Three different variants of such methods 
have been developed that combine the results from the 
different short segments incoherently: the "stack-slide", 
the "Hough-transform" and the "PowerFlux" schemes. 
The stack-slide procedure jTSJ [5U] averages the normal- 
ized power from the Fourier transform of 30-minute seg- 
ments (Short time baseline Fourier Transform, SFT for 
short) of the calibrated detector strain data. The Pow- 
erFlux scheme [53[ [25] can be seen as a variation of the 
stack-slide method, where the power is weighted before 
summing. The weights are chosen according to the detec- 
tor noise and antenna pattern to maximize the signal-to- 
noise ratio (SNR). The Hough-transform method [1^1127] 
sums weighted binary counts, depending upon whether 
the normalized power in an SFT bin exceeds a certain 
threshold. 

As the segment duration is increased, it becomes nec- 
essary to account for signal modulations within each seg- 
ment by computing the ^-statistic over a grid in the 
space of phase evolution parameters. This results in a sig- 
nificant increase in the computational requirements of the 
search. The distributed volunteer computing project Ein- 
stein@Home [28] has been created to address this need. 
Two previous papers [23] [21] report on results of such 
CW searches from the fourth LIGO science run and from 
the first two months of S5, respectively. The method used 
was based on the computation of the coherent ^-statistic 
on data segments from either the H or L detectors sepa- 
rately, and only parameter space points with values of 2F 
larger than 25 were returned back to the Einstein@Home 
server for further inspection. The threshold value of 25 
limited the ultimate sensitivity of that search: if a sig- 
nal was not loud enough to surpass that threshold on at 
least some of the segments it would not be detected. The 
threshold value was set by bandwidth constraints on the 
size of the results file uploaded back to the server by the 



6 



host, i.e. on the maximum number of significant points 
that could be returned. These results were subsequently 
combined by a coincidence scheme, performed offline in 
the post-processing phase. In contrast, in the search pre- 
sented here, the combination of the results from the co- 
herent searches takes place directly on the host machines 
using a Hough-transform scheme. This makes it possible 
to use a much lower threshold on 2T, equal to 5.2, that 
defines the parameter space points to be passed on to the 
Hough-transform. Moreover, in this search, data from 
the H and L detectors are coherently combined O [H] . 
Finally, more data was searched in this analysis com- 
pared to any previous Einstein@Home search. The Ein- 
stein@Home runs presented here refer to searches based 
on the first (S5R3) and second year of S5 LIGO data. 
This latter search was run on Einstein@Home in two sep- 
arate steps, called S5R5 and S5R6. Since the S5R6 run 
used the same data as S5R5, but extended the search re- 
gion above 1 kHz, in this paper we simply refer to these 
two runs as S5R5. 

The paper is structured as follows. In Sec. |Hj we 
briefly review the Hough-transform method. Section |III| 
describes the Einstein@Home distributed search used to 
analyze the data set. Section [TV] gives a detailed descrip- 
tion of the S5R5 post-processing, which is based on the 
pioneer S5R3 post-processing (described in Appendix[B|). 
Upper limit computations from the more-sensitive S5R5 
data are provided in Sec.[V] The study of some hardware- 
injected signals is presented in Sec. |VI| In Sec. |VII| we 
make some concluding remarks. 



II. THE DATA ANALYSIS METHOD 

A. The waveform model 

Let us begin by briefly describing the standard signal 
model for CW signals. In the rest frame of the neutron 
star, the gravitational wave signal is elliptically polarized 
with constant amplitudes A+. x for the two polarizations 
h + x (t). Thus, we can find a frame such that 

h+(t) = A + cos<p(t) , h x (t) = A x sin</>(t). (1) 

The two amplitudes are related to an overall amplitude 
ho and the inclination angle t between the line of sight 
to the neutron star and its rotation axis 
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+ = -/io(l + cos l) , A x =/i cost. (2) 



The value of h is model dependent. For emission due to 
non-axisymmetric distortions, the amplitude ho depends 
on the ellipticity e of the star defined as 



\Ixx Iyy\ 



(3) 



Here I zz is the principal moment of inertia of the star, 
and I xx and I yy are the moments of inertia about the 



other axes. The amplitude is given by 

_4ir 2 GI zz f 2 e 
n - — -a ; — , 



(4) 



where / is the frequency of the emitted GW signal (which 
is also twice the rotational frequency of the star), G is 
Newton's constant, c is the speed of light and d is the 
distance to the star. The distribution of e for neutron 
stars is uncertain and model dependent since the break- 
ing strain for a neutron star crust is highly uncertain (see 
e.g. [21 I29H3T] for further discussion). 

Energy loss from the emission of gravitational and/or 
electromagnetic waves, as well as possible local acceler- 
ation of the source, causes the signal frequency arriving 
at the solar system to evolve. To first order, it can be 
expressed as 



f{r) =k + f(r ~ T ), 



(5) 



where r is the arrival time of a wavefront at the solar sys- 
tem barycenter (SSB), /o is the frequency at a fiducial 
reference time To, and / denotes the first time derivative 
of the frequency. The astrophysical implications of ig- 
noring higher order derivatives in this Taylor expansion 
will be discussed later. The phase of the signal, <j>(t), fol- 
lows directly from the frequency evolution with an initial 
phase 4>o at the reference time. 

As the detector on the Earth moves relative to the SSB, 
the arrival time of a wavefront at the detector, t, differs 
from the SSB time t: 

f{t) ■ ft 



r(t) =t + 



EQ 



(6) 



Here f(i) is the position vector of the detector in the 
SSB frame, and ft is the unit vector pointing to the neu- 
tron star; Aeq and Ag are respectively the relativis- 
tic Einstein and Shapiro time delays [35]. In standard 
equatorial coordinates with right ascension a and decli- 
nation 5, the components of the unit vector ft are given 
by (cos a cos 5, sin a cos 5, sin<5). 

Ignoring the relativistic corrections, the instantaneous 
frequency f(t) of a CW signal, as observed at time t 
by a detector on Earth, is described by the well-known 
Dopplcr shift equation: 



/(*) = />) + />) 



v{t) 



(7) 



where v(t) is the detector velocity with respect to the 
SSB frame; v(t) is the sum of two components, from the 
yearly Earth motion around the Sun (v y ) and from the 
rotation of Earth around its axis (vd)- 

Finally, the received signal at the detector is 

h(t) = F + {t;ft,iP)h+{t) + F x {t;ft,iP)h x (t), (8) 

where F + x are the detector beam pattern functions 
which depend on the sky position ft and the relative po- 
larization angle ip of the wave-frame (T2l [33] . There are 
thus altogether eight signal parameters, which include 
the four phase evolution parameters (/o, /, a, 5), and four 
other parameters (ho,L,tp,<j)o). 
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B. The Hough-transform algorithm 

For completeness, we summarize the Hough detection 
statistic in this section. Further details of the method 
are given in [Tj5] and previous searches with this method 
applied to short coherent times are reported in [2TJ [22] ■ 

The Hough-transform is a well-known technique used 
mainly in digital image processing for robust extraction 
of patterns. Such a procedure is employed here for iden- 
tifying points in the time-frequency plane that match the 
pattern expected from a signal. The time-frequency data 
in our case is the /-statistic computed as a function of 
signal frequency, for each of the data segments of dura- 
tion T seg , over a grid of points in the space of (a,S,f). 
The grid in (a,8,f) space used for this /"-statistic com- 
putation is called the coarse grid because its resolution is 
determined by the coherent time baseline T seg , and makes 
no reference to the full observation time or the number 
of segments N scg . The result of this computation is thus 
a collection of /"-statistic values T % . -(k) where the in- 

tegers k and i label a frequency bin (with spacing Sf as 
defined below) and a data segment respectively. 

The frequency and frequency derivative spacings for 
the coarse grid are based on choosing the maximum al- 
lowed fractional loss in the /-statistic when the signal 
and template points are slightly mismatched. This leads 
naturally to the notion of a metric in parameter space 
[5H[3S] and this has been studied for the CW case in [13] . 
The grid spacings in /, / are given respectively by [171 131)] 



Sf = 



y/l2r 



and 



Sf 



ttT 2 

scg 



(9) 



(10) 



where m represents the nominal single dimension mis- 
match value. For all the Einstein@Home runs discussed 
here m has been taken equal to 0.3. The span T seg 
of each segment has been set equal to 25 hours for all 
the runs. The frequency resolution, given by Eq. Q, is 
Sf ~ 6.7 /iHz for all the Einstein@Home runs described 
in this paper. As we shall see shortly, for technical rea- 
sons it turned out to be necessary to use a finer spacing 
for / than given by Eq. (10). 

In combining these N se g different /-statistic vectors, 
it is necessary to use a finer grid in (a, 6, f) centered 
around each coarse grid point. Our implementation of 
the Hough-transform algorithm assumes that the fine 
grid is a Cartesian product of a rectangular sky-grid and 
one dimensional grids in / and /. Moreover the fine sky- 
grid is assumed to be aligned with the a and S directions 
in the sky. In order to completely cover the sky with 
the different fine sky-grids, it is thus simplest to choose 
a rectangular coarse sky-grid aligned with the (a, 5) di- 
rections. We choose a coarse grid such that the spacing 



in 6 is a constant and the spacing in a is proportional to 
(cos<5) -1 . This ensures that each cell of the coarse grid 
covers a fixed solid angle. 

Since the coherent integration time is very close to a 
sidereal day, it is reasonable to assume that the Hough 
sky-patch size d6 (which in our case is the same as 
the coarse sky-grid resolution) is determined by Vd, the 
Earth's rotation speed at the equator. At frequency / 
we have H91: 



de 



(11) 



.scg 



In practice, this estimate was verified by Monte-Carlo 
studies with signal injections, and we used 



dO 



d6 for S5R3 
Tide for S5R5, 



(12) 



where the factor 7Z = y3 increases the size of the sky- 
patch d9 for the S5R5 run. The Monte-Carlo studies 
showed that the grid spacing for S5R5 corresponded ap- 
proximately to a mismatch of m ~ 0.3 so that, for any 
other value of m, the spacing would be approximately 
■s/m/O.ZKde. 

We also need to set the resolution of the refined sky- 
grid used by the Hough algorithm. Since the full obser- 
vation time is of the order of a year, the relevant scale 
here is set by the speed v y of Earth as it orbits the Sun. 
Following [19], the resolution for the fine sky-grid at a 
frequency / is given by 



cSf 

PfVy 



(13) 



The parameter p scales the resolution compared to the 
conservative estimate cSf/(fv y ) and in practice, again 
based on Monte-Carlo studies, we used p = 0.5. One can 
see that the increase in the number of sky position points 
from the coarse /-statistic grid to the fine Hough grid is 



ref 
sky 



(d^/dfln) 



\VV2m v d 



~O(10 4 ). (14) 



Taking m = 0.3 and p = 0.5 yields Af^ ~ 8444 for the 
S5R5 run. 

Finally, let us turn to the coarse and fine grids for 
/. Ideally, we should refine the coarse / grid spacing of 
Eq. (10) by a factor N seg [T5]. However, in our imple- 



mentation of the Hough transform, using this refinement 
turns out to increase the maximum memory footprint of 
the different searches and would make it unsuitable for 
EinsteinQHome. As a compromise, it was decided not to 
use any refinement in / and instead to use a finer resolu- 
tion for the coarse grid. Based on Monte-Carlo analyses, 
an acceptable value for the / spacing turns out to be 



in 



rp2 



1.2 X 10 



-10 



Hz 



(15) 



This value was used for S5R5. However, for S5R3, this 
was incorrectly set to 



(Sf 



S r .l;:i 



J^2 



3.8 x 1CT 10 Hz s 



(16) 



leading to a corresponding loss in sensitivity for S5R3. 

The flow chart of the search algorithm used for this 
search is depicted in Fig. [T] The input data set is com- 
posed of 30-minute baseline SFTs. This set is partitioned 
in subsets such that no more than 25 hours of data are 
spanned by each segment and such that there is overall 
(including data from both detectors) at least 40 hours of 
data in each segment. Let T D b s be the observation time 
spanned by the N scg segments constructed in this wajQ 
The multi-detector .F-statistic is computed for each seg- 
ment at each point of the search parameter space (/,/, a, 
S). The next step consists of selecting parameter space 
points for which the /"-statistic is above the fixed thresh- 
old. For every set of (/,/, a, 5), we assign a value n; = 1 
or in the i th segment depending on whether the cor- 
responding .F-statistic is above the threshold 5.2 or not; 
this threshold turns out to be optimal [TJJ]. The values 
n-i(f) are called a "peakgram", which is the input to the 
Hough-transform. 

The final statistic used by the Hough search is a 
weighted sum of binary counts Uj, giving the so-called 
Hough number count n c , expressed by [22j 



JV B „ g -l 

E < 

i=0 



(17) 



The weight Wi for a frequency / and for a particular sky- 
location is determined from the average antenna response 
and average detector noise over the duration of the i th 
segment. Since the input data in each segment consists 
of SFTs, we perform the averaging over each SFT. Let 
Ni be the number of SFTs in the i th segment. Let S]^ 1 
be the single-sided power spectral density (PSD) of the 
7 th SFT in the i th segment, and averaged over a narrow 
frequency band containing the search frequency. The Wi 
are given by 



Ni-l 



Wi (X 



E (*+«.,) 



x(vy) J Wi >"> ' 



(18) 



7=0 



where F+u^) and F x u^\ are the detector antenna pat- 
tern functions for the 7 th SFT in the i th segment, and 



'ft \ 1 0=0 b h J 



(19) 



Produce segments of data spanning 
no more than 25 hours and containing 
at least 40 hours of total data 



^Compute the F-statistic 
for each segment 



For everv [txjftA ) select 
frequency bins and 
produce peakgrams 



Perform the Hough _^ Return top 10 000 .1 

~* j.j . .,„ , — H Off-line post-processing 
transform candidates per WU I - 



FIG. 1. High level schematic of the pipeline used in 
the searches. For each data segment, the multi-detector T- 
statistic is computed and frequency bins are selected setting 
a threshold on the .F-statistic. Such selected frequency bins 
are then used to create the Hough map. The output is a 
set of candidates in the parameter space. The color-filled 
boxes indicate the steps performed on the volunteer comput- 
ers. The acronym "WU" is defined in Sec. |III A] and refers to 
the independent computing tasks into which we partition the 
computational work of the search. 



It is easy to see that, in the hypothetical case when the 
data is exactly stationary, so all the S^f 1 are identically 
equal to each other, then u^ j7 = 1. More realistically, the 
u>i n are approximately unity for stationary data and the 



use of the harmonic mean in Eq. ( 19 ) ensures that the 



z j7 



do not deviate too far from unity in the presence of 



non-stationary noise. 

The weight normalization is 



Wscg-l 

E < 

i=0 



scg , 



(20) 



which ensures that the Hough number count n c lies 
within the range [0, jV seg ]. In [22] it is shown that the 
weights Wi, first derived in |37| . maximize the sensitivity, 
averaged over the orientation of the source (see Appendix 
|A| for a further discussion of the weights and some tech- 
nical problems that were encountered in the search). 

From the Hough number count n c we define the signif- 
icance (or critical ratio), CR: 



CR 



(21) 



that measures the significance of the measured n c value 
as the deviation from the expected value n c = N seg p 
in absence of any signal, in units of the expected noise 
fluctuations a; pis the probability that a parameter space 
pixel is selected in the absence of a signal. In case of unity 
weighting, the standard deviation is simply that of the 
binomial distribution: a = U N scg p(l — p). When the 
weights are used, the standard deviation is given by 



1 Note that T seg ^ Tobs/TVseg because of the gaps that are un- 
avoidably present in the data stream, corresponding to times 
when the interferometers were not in lock, and to the selection 
of 25 hour segments containing the requisite amount of data. 



o- = y/\\w\\*p(l-p), 



(22) 



where 



= E 



JVs, 



-1 2 

wf 



The CR is the detection 



statistic returned by the hierarchical searches presented 
here. 
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III. THE EINSTEINQHOME DISTRIBUTED 
SEARCH 

The Einstein@Home project is built upon the BOINC 
(Berkeley Open Infrastructure for Network Computing) 
architecture [3&H40] , a system that exploits the idle time 
on volunteer computers to solve scientific problems that 
require large amounts of computer power. During the 
S5R5 run, Einstein@Home had approximately 225 000 
registered volunteers and approximately 750 000 regis- 
tered host machines that contributed a total of approxi- 
mately 25 000 CPU (Central Processing Unit) years. 



A. Details of the S5R5 search 

The computational load is partitioned in independent 
computing tasks, called "workunits" (WUs), each of 
which is analyzed by a volunteer machine. In particu- 
lar, 7 369 434 and 10 945 955 WUs have been generated 
for the S5R3 and S5R5 runs, respectively. The S5R5 
run was launched on January 13, 2009 and ended on 
October 30, 2009. It used 10 560 S5 LIGO SFTs, col- 
lected between the GPS times of 852 443 819 s (Wed Jan 
10 05:56:45 GMT 2007) and 875 278 812 s (Mon Oct 01 
12:59:58 GMT 2007). The analyzed data consists of 5 550 
and 5 010 SFTs from the LIGO H and L interferome- 
ters, respectively. The number of data segments used 
for the S5R5 run is 121, each spanning no more than 
T scg = 25 hours and with at least 40 hours of data, as al- 
ready said. Similar details for the S5R3 run can be found 
in Appendix |B| 

The total search frequency range of the S5R5 search is 
[50, 1 190] Hz, with a frequency resolution Sf ~ 6.7^iHz 
and spindown resolution Sf ~ 0.12 nHz s _1 . Each WU 
analyzes a constant frequency band B ~ 20 mHz, the full 
spindown interval, ranging roughly from —2 nHz s _1 to 
0.11 nHz s _1 , and a region of the sky, as we shall see in 
Sec. lHTCl 

The original data contained instrumental artifacts in 
narrow frequency bands that were known before the 
launch of the Einstein@Home run. Those bands were 
identified and the corresponding frequency bins in the 
SFTs were replaced with white Gaussian noise at the 
same level as the neighboring frequencies. Table [I] shows 
which bands were treated in this manner and what in- 
strumental artifact they harbored. These control bands 
are useful to compare and contrast the results obtained 
on real data against pure theoretical noise. Measure- 
ments and studies after the Einstein@Home run refined 
the frequencies and widths of these artifacts and identi- 
fied additional ones; the final lists of artifacts are given in 
Appendix [Cj and were used to discard candidates (as we 
shall see in Sec. IV B). The "cleaning" process affected 
ih bandwidth, in addition to bands that 



Cause 



/l (Hz) Harmonics LFS (Hz) HFS (Hz) IFO 



IVBJ). The 

~ 27 Hz of searcl 
were eliminated later in post-processing. 

The output data files from each WU are stored as ZIP- 
compressed ASCII text files containing the 10 000 most 
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TABLE I. Instrumental lines identified and "cleaned" be- 
fore the EinsteinQHome runs. The different columns repre- 
sent: (I) the source of the line; (II) the central frequency of 
the instrumental line; (III) the number of harmonics; (IV) 
Low-Frequency-Side (LFS) of the knockout band; (V) High- 
Frequency-Side (HFS) of the knockout band; (VI) the inter- 
ferometer where the instrumental lines were identified. Note 
that when there are higher harmonics, the knockout band 
width remains constant. 



significant candidates ranked according to the signifi- 
cance (as defined in Eq. (21)) over the parameter space 
searched by that WU. The decision to keep the top 10 000 
candidates was based on the maximum upload volume 
from the hosts to the Einstein@Home servers. All in all, 
on the order of 10 11 candidates were returned to the Ein- 
stein@Home server from each run, corresponding roughly 
to 2.3 TB of data. 



The files contain nine quantities for each candidate. 
The first four are the values (on the coarse grids) of the 
frequency /, sky-position (a, 5), and spindown /. The 
fifth is the significance of the candidate as defined in 
Eq. (21). The remaining four quantities are connected 



with the refined sky-grid centered on each coarse sky- 
grid point (recall that refinement is performed only on 
the sky). In particular, the location of the most signifi- 
cant point on the fine sky-grid, and the mean and stan- 
dard deviation of the Hough number count values on all 
points of the fine grid are returned. 
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B. Validation of returned candidates 

In order to eliminate potential errors, due to defective 
hardware and/or software or to fraud, BOINC is con- 
figured so that each WU is processed redundantly by 
computers owned by at least two different volunteers. 
An automated validation process checks the consistency 
of the results, ruling out those that are inconsistent, in 
which case new WUs are generated to run again indepen- 
dently. The first step of the validation is to check that 
the file syntax is correct and that the first four values, 
i.e. (/, a, S, /), are within the appropriate ranges. Next, 
for the pair of result files from each WU, the validator 
checks that the values of (/, a, S, f) agree to within float- 
ing point accuracy (the frequency is in fact checked to 
double precision). Finally, the significance values CRi 
and CR2 from the two result files are compared and are 
validated if 



|CRi - CR 2 |/(CRi + CR 2 ) < 0.12 . 



(23) 



In S5R5, about 0.045% of results that were processed by 
the validator were marked as "invalid", including both 
syntax errors in individual files and errors in comparisons 
of different results files. Excluding the syntax errors in 
individual files, the error rate arising from comparisons of 
pairs of distinct result files (most likely due to differences 
in floating point arithmetic on different computational 
platforms) was ~ 0.015%. 

Is it possible that two invalid results could agree with 
each other and thus end up being marked as valid? While 
it is difficult to exclude this scenario with complete cer- 
tainty, an upper limit for the probability of this happen- 
ing is (0.015/100) 2 w 2.2 x 10~ 8 (only the 0.015% error 
rate due to comparisons of distinct result files is rele- 
vant here). As mentioned earlier, there were a total of 
~ 1.1 x 10 7 WUs. It is therefore unlikely that even a 
single pair of result files would be incorrect and still pass 
validation. 

The threshold of 0.12 on the value of A defined above 
turns out to be much looser than necessary. The differ- 
ences in the actual observed values of A from a pair of 
matching result files are usually much smaller. The ob- 
served standard deviation of A turns out to be ~ 0.012, 
i.e. an order of magnitude smaller than the threshold. 
In addition, the standard deviation of the difference 
CRi — CR 2 is measured to be ~ 0.15, which corresponds 
to a standard deviation of ~ 0.7 in terms of the number 
count. As we shall see later (see e.g. Fig. [3]), the loudest 
events in every 0.5 Hz band have an average loudest num- 
ber count > 70. Thus, these differences correspond to a 
< 1% effect in the number count at the 1-cr level, and we 
expect this to have a negligible effect on our analysis. 



each WU must cover the same number of parameter space 
points so that they can be completed in roughly the same 
amount of time by a typical host machine. Second, the 
amount of data that must be downloaded by each host 
machine and the maximum memory footprint of each job 
must be within appropriate limits. Finally, one needs to 
choose the computational time for each WU on a typical 
host machine and the total time that the project should 
run, given the total computational power that is avail- 
able. To meet these requirements, we need to understand 
how to split up the parameter space, and to measure the 
CPU core time spent by the search code on each part of 
the analysis. 

We start with the basic parameters of the search, 
namely the total observation time T b s , the coherent time 
baseline T seg and the number of segments N seg , the res- 
olution for the coarse and fine sky-grids, dOjr, d9n, given 
by Eqs. ( 12 1 and ( 13 ), respectively and the frequency and 
spindown resolutions, Sf and Sf, given by Eqs. ^ and 



( 15 ), respectively. Recall that the limit on the maximum 



memory footprint of each job already forced us to forego 
any refinement in /. 

Unlike the previous Einstein@Home search [Mj , where 
each WU searched the whole sky, here we choose each 
WU to cover a fixed frequency bandwidth B, the entire 
spindown search range, and a limited area of the sky. Let 
fb be the highest frequency in the 6 th search band. We 
want the computation time for every WU to be approxi- 
mately the same, hence every WU must search the same 
number of coarse sky-grid points. Since the resolution in 
the sky, dOjr, is inversely proportional to the frequency 
of the signal that we are searching for, WUs at higher 
frequencies will be searching smaller portions of the sky. 
Let Nb be the number of coarse grid sky points over the 
whole sky for the 6 th band: 



47T 



(?) 



2 f, 2 T 2 

Jo scg 



(24) 



which is shown in Fig. [2] for the S5R5 run. At frequency 
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C. Workunit design 

The design of the WUs must satisfy certain require- 
ments. The first is that the WUs must be balanced, i.e. 



FIG. 2. Total number of J-"-sky-grid points, for the S5R5 
run, as a function of the frequency, given by Eq. (24 1. 



fb, we then partition the sky in Pb parts, each containing 
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iVsky points: 



with K = Air i 



Ph 



K 



sky 



(25) 



' T*JRr*K 



sky 



In practice, to limit 
the number of sky-grid files needing to be downloaded 
by the host machines, the sky-grids are constant over 
10-Hz frequency bands and are determined based on the 
highest frequency in each band. 

We choose a constant range of / values over the entire 
frequency band. To see what this implies for potential 
sources, we define the spindown age for a system emitting 
at a frequency / and spindown / to be r = f /\f\ [^J It 
is clear that at a given search frequency, the minimum 
spindown age is determined by the maximum value of |/| 
included in the search. Our choice of a constant range 
of / means that the minimum spindown age is frequency 
dependent: 



h 



l/lr 



(26) 



We choose the minimum age at 50 Hz to be 800 years, 
and a fixed range in /: ~ [—20, 1.1] x 10~ 10 Hz s _1 . 

The total number of WUs is simply the sum of the 
total number of sky-partitions at each frequency //,: 



N 



wu 



K 
~B 

K 
3B 



f 2 df 



(/max /min) • (27) 



In the third step we have replaced the sum over frequen- 
cies with an integral between a minimum (/ m i n ) and max- 
imum (/ ma x) frequency. As a consistency check, we shall 
see later that B ~ 20mHz, which is sufficiently small 
that this is a good approximation. 

The total computing time for one WU can be expressed 

as 

rwu = NtkyNf [{N f + N sh )T? + N f T H ] , (28) 

where Nf , Nf , and iV s k y represent the number of coarse 
search frequency bins, spindown values, and sky-grid- 
points respectively, Tjr and th are the times needed to 
compute the ^-statistics and Hough number count for 
one point of the coarse grid; N s \, ^\ represents addi- 
tional "sideband" bins needed to compute the Hough- 
transform. The need for these sidebands can be under- 
stood by thinking about computing the Hough number 



count for a frequency near the edge of a given search 
band: that number count will involve summing the peak- 
grams along a curved track that can extend a small dis- 
tance to cither side of the target frequency. 

The time Tjr needed to compute the ^-statistic for one 
point of the coarse grid can be expressed as 



,jrN S 



FT, 



(29) 



where N$ft is the total number of SFTs used in the 
search; Tjr is the time needed to compute the ^-statistic 
per single SFT. The time th in Eq. (28), needed to com- 



pute the Hough number count corresponding to one point 
of the coarse grid, can be written as 



H-^seg-A/'sky, 



(30) 



where is the time to sum the Hough number count 
per single data segment and per single point of the fine 
grid. Here is an overall refinement factor, i.e. the 

number of grid points analyzed by the Hough algorithm 
for each coarse grid point. In our case, since the only 



refinement is over the sky, Af^ y is given by Eq. (14). 
The computational time is thus determined by the two 
timing constants Tjr and r^. For our implementation of 
the algorithm, these constants were measured to be 



Tj- = 180 ns, t„ = 1.1 ns 



H 



(31) 



These numbers are of course only average values for a 
typical host CPU core available at the time of the Ein- 
stein@Home runs. 

The presence of N s ^ leads to an overhead for the com- 
putation. We want to control this overhead and keep it 
below some acceptable level. Thus, we define the over- 
head e to be the ratio between the time spent in com- 
puting the ^-statistic for the "sidebands" and the total 
computational time: 



A/ S b tj- 



(32) 



The bandwidth B needs to be sufficiently large so that 
e is sufficiently small, but a too large value of B can 
lead to excessively high download volumes for the Ein- 
stein@Home clients. From the above equation, we see 
that the frequency bandwidth B can be determined by 
fixing e: 



B = 



A/ sb Tjr fl-e 



(33) 



2 The reader is warned that the spindown age as defined here can 
be rather different from the true chronological age of the star. 
For a neutron star with a present day frequency / and spindown 
/ due purely to the emission of GWs, the age would be //4|/|; 
more generally this depends on the physical mechanisms that are 
responsible for the spindown |41) . 

3 The bins 7V sb can be calculated by A/ ab /<5/, where the average 
Hough "sidebands" (A/ S t,) can be estimated from the Hough 



master equation. The frequency offset corresponding to the half 
diagonal distance An, over one Hough sky-patch, is (A/ S b) = 

t{\v y ■ An\) = ■ By using T scg = 25 hours, we get for 

S5R3 (A/ sb ) ~ 5 x 10~ 4 Hz and N Bb ~ 75 bins, while for S5R5 
(A/ sb > ~ 8.7 X 1CT 4 Hz and 7V sb ~ 130 bins. However, to be 
conservative, (2 X iV s b) bins were used for safety reasons and to 
take into account for the changes of the velocity on the different 
segments. 



12 



For all the Einstein@Home runs presented here we choose 
e = 5%. For the S5R5 run, this leads to B ~ 20 mHz. 
The total run-time r p of the project is 



T"WU %U 



N, 



CPU 



(34) 



where -ZVcpu represents the number of volunteer CPU 
cores. Given a certain number of CPU cores and hav- 
ing fixed T p , the maximum search frequency / max can be 
derived from the above equations to be 



3 

max 



J mm 1 



^TpNcpu (1 



e) 6f 



K N 



(35) 



/ 



where K = 4tt (^) T s 2 eg ^~ 2 . For the S5R3 and S5R5 
runs, the nominal project duration was chosen to be 6 
months. With the above choices, the search frequency 
ranges for S5R3 and S5R5 turn out to be respectively 
[50, 1 200] Hz and [50, 1 190] Hz. 



D. Accuracy of spindown model 

Let us briefly discuss the second order spindown / 
which, as mentioned previously, is not a part of our 
search. For our frequency resolution 8f, given by Eq. ^, 
and the full observation time T b s , / would have to be at 
least 



fn 



8f/T< 



obs 



1.3 x 10" 2O Hz s 



(36) 



in order for the signal to move by a single frequency bin 
over the full observation time. On the other hand, for 
a minimum spindown age T m i n at a frequency /, a use- 
ful estimate for the range of / that we should search is 
// T min = /max//- Th us, / is potentially more important 
at lower frequencies and higher spindown values. Using 
the maximum value of |/| in the search and the minimum 
value of the search frequency gives us a value of / that 
might be potentially of astrophysical interest: 



/ast 



x 10 _20 Hz s -2 



l/l* 



2.0 x 10~ 9 Hz/s 



50 Hz 



(37) 



Comparing with the minimum value of / obtained above, 
we see that there is potentially a region in (/, /) space 
where we could improve our astrophysical detection ef- 
ficiency by including /; for our chosen range of / there 
is no effect of / above ~308 Hz. It is important to note 



that the calculation of Eq. (36) is too conservative be- 



cause it does not include any correlations between the 
phase evolution parameters. On the other hand, there 
is considerable uncertainty in the value of / as t- If the 
neutron star has a braking index n = ff/f 2 , then / ast 
increases by a factor n. If a star is spinning down purely 
due to gravitational wave emission, then n = 5. On the 



other hand, for the Vela and Crab pulsars, observed val- 
ues of n are ~ 1.4 and ~ 2.5 respectively [42| 143] . 

The actual impact on our astrophysical reach is thus 
hard to quantify. Let us consider as an example the ex- 
treme case when the spindown is entirely due to grav- 
itational radiation so that the braking index is n = 5. 
For this case, a conservative estimate of the part of the 
spindown range |/| cons that is included in our search is: 



2 x 10" 9 Hz s" 1 



< 0.18 



50 Hz 



(38) 



Here we have used Eqs. (36) and (37) modified by the 
braking index factor, 
spindown value of 
of 1190 Hz and 



This corresponds to a minimum 
T.8 nHz s _1 at the upper frequency 
-3.6 x 10- 10 Hz s- 1 at 50 Hz. 



IV. 



S5R5 POST-PROCESSING 



As said earlier, roughly 10 11 candidates from the S5R5 
run were returned to the Einstcin@Home server. They 
were then transferred to the 6720-CPU-core Atlas Com- 
puting Cluster [44 at the Albert Einstein Institute in 
Hannover, and post-processed. The goal is to filter the 
set of 10 11 candidates, excluding false candidate events. 
The post-processing strategy consists of the following 
steps: 

• selection of 100 most significant candidates in 
0.5 Hz frequency bands; 

• removal of known instrumental noise artifacts; 

• removal of unknown data artifacts through the T- 
statistic consistency veto; 

• follow-up of the most significant candidates with 
S5R3 data; 

• fully-coherent follow-up of the surviving candi- 
dates. 

The items outlined above are described in the next sub- 
sections. 



A. Selecting the top candidates in frequency bands 

As is commonly done in CW searches, the results are 
examined separately in fixed-size search frequency bands; 
here we choose to perform the analysis in 0.5 Hz bands. 
As described earlier, in designing the WUs, we have 
previously been led to break up the frequency range in 
~ 20 mHz bands and the sky has been partitioned as well. 
This was however done for purely technical reasons to 
make the search on Einstein@Home feasible. The choice 
of frequency bands for the post-processing is based on dif- 
ferent requirements. First, we would like the detector to 
have roughly constant sensitivity within each frequency 
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band. Furthermore, as we shall see, the search does not 
result in a convincing detection candidate, and upper lim- 
its will be set over each of these frequency bands. Hav- 
ing a large number of very narrow bands would make the 
calculation of the upper limit very computationally in- 
tensive. The choice of 0.5 Hz is a compromise between 
these two requirements. This choice is in fact compara- 
ble to previous CW searches and will make comparisons 
straightforward. Finally, we note that all other things 
being equal, having a larger frequency band will in prin- 
ciple also lead to a decrease in sensitivity simply because 
of having a larger number of templates. This is however 
a relatively minor effect in the present case. 

For each of the 0.5 Hz- wide frequency bands, we select 
the 100 most significant candidates for further analysis, 
leading to a set of 228 000 loudest S5R5 candidates. This 
choice was dictated by the available computational and 
human resources for the post-processing. As will be il- 
lustrated in the following, at the end of the automated 
post-processing procedure there will remain of order 10 
candidates which survive all selection criteria. This is 
about the number that we can afford to follow-up man- 
ually with further investigations. As our follow-up pro- 
cedures are further automated and optimized, it will be- 
come possible to consider lower thresholds and to inspect 
a correspondingly larger number of candidates. 

The number count of such candidates is plotted in 
Fig.[3]as a function of the frequency. For most bands, it is 
consistent with expectations. The number count gener- 
ally increases with increasing frequency because the num- 
ber of sky points searched over increases (see Eq. (11)) 



and the maximum expected value of a random variable, 
over repeated trials, grows with the number of indepen- 
dent trials. Let us ignore the effect of the weights, and 
take the Hough number count n c (defined in Eq. (17)) to 
be an integer random variable following a binomial dis- 
tribution. The cumulative probability for obtaining n c 
or lower is 



F(n c ) = 



scg 

n 



(39) 



The binomial parameters are N scg = 121 and p scg = 
0.267 (consistent with a 2J- threshold of 5.2). The prob- 
ability p max that the maximum number count is n c over 
a set of At r i a i s independent trials is 



Atrial, 



Pmax(n c , NtrialslNseg^seg) = F{n c ) Nt ^ ls -F(n c -1) 

(40) 

We compute Eq. (40) as a function of N tl i a i s , which we 



take to be the number of templates searched to cover 0.5 
Hz bands, i.e. 7.5 x 10 4 frequency values x 18 spindown 
values x 8 444 Hough pixels per T sky point x Nb, the to- 
tal number of T sky points shown in Fig. [2j The number 
of Hough pixels has been computed using Eq. ( 14 ) with 
m = 0.3 and p = 0.5. Figure [3] shows the expected value 
of the maximum (central black curve ), c omputed using 
the probability function given by Eq. ( 40 ) , superimposed 



there is broad agreement, in most bands, between our 
results and the expectations for Gaussian noise. 



B. Removing known data artifacts 

As a first step of the post-processing pipeline, we elimi- 
nated from the list of top candidates any candidate whose 
frequency was too close to that of either a known artifact 
or to the cleaned noise bands described in section MI Al 
Specifically, we discarded those candidates whose detec- 
tion statistic could have been constructed with contribu- 
tions either from: 

• data polluted in either of the two instruments by 
spurious disturbances; details of such detector dis- 
turbances are given in Appendix [C] and, in particu- 
lar, a list of known spectral disturbances for the H 
and L instruments are listed in Tables I VI I and IVHI 

• from fake noise that had been inserted by the clean- 
ing process and hence could not host a CW signal 
(see Table [T] in Sec.|lHA]). 



on our measurements (red circles) and confirming that 



After this veto, about 25% of the candidates were elimi- 
nated from the original set of 228 000 loudest S5R5 can- 
didates. More precisely, a total of 172 038 S5R5 candi- 
dates survived this veto. The bandwidth removed due to 
the lines listed in Appendix [C] amounted to ^243 Hz; an 
additional 27 Hz was removed due to the cleaned noise 
bands. 



C. The ^-statistic consistency veto 

We have thus far considered only known instrumental 
disturbances for vetoing candidates. However, we expect 
there to be more such disturbances present in the data 
that have not yet been explicitly identified. The idea is 
to discriminate between disturbances in a single detector 
and signals, which should produce consistent values of 
the ^-statistic in both detectors [45] . We refer to this 
method as the ^-statistic consistency veto. 

For each of the 172 038 S5R5 surviving candidates, the 
single-detector and multi-detector ^"-statistic was com- 
puted for each of the 121 data segments and then aver- 
aged over the segments. We refer to these averaged 2T 
values as (27^) and {2F\) for the H and L detectors 
respectively, and (2J 7 hl) for the coherent combination 
of the data from the two detectors. Candidates were dis- 
carded if either (2.Fh) or (2Jx) were greater than (2J 7 hl)- 
Using this veto, a small fraction (4.1%) of candidates was 
eliminated, leaving 164 971 surviving S5R5 candidates. 

The impact of this veto is limited in this case due to the 
prior removal of the bulk of instrumental artifacts. How- 
ever, the ^-statistic consistency veto represents an effi- 
cient method to remove disturbances that clearly stand 
out of the noise in the absence of independent instru- 
mental evidence. Figure [4] shows the average 2J r -values 
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FIG. 3. Number count of top 100 loudest candidates (blue dots) selected in 0.5 Hz-wide frequency bands as a function of the 
search frequency, across the entire S5R5 search frequency range. The loudest (most significant) candidate in every 0.5 Hz band 
is indicated by a red circle. The expected values of the loudest candidates for Gaussian noise alone are shown by the central 
black curve. The lower and upper black curves show ±3 standard deviations from the expected value. 



for 164 971 surviving (top plot) and 7 067 vetoed (bottom 
plot) S5R5 candidates as a function of the multi-detector 
average 2J r -values. By construction, all the surviving 
candidates in the top panel of Fig. [4] lie below the red 
dotted line, which defines the veto criterion. 



D. Distribution of candidates 

We have now applied all of our vetoes that try 
to remove instrumental artifacts. While there will of 
course remain other low amplitude instrumental spectral 
lines and hardware signal injections (described later in 
Sec. VI), we now need to deal with the possibility that, 



say, even Gaussian noise can mimic a signal in some cases. 
All remaining candidates will need to undergo detailed 
individual inspection and we will only be able to afford 
this for a few candidates. As our follow-up techniques be- 
come more refined, optimized and automated, we will be 



able to improve this part of the pipeline and dig deeper 
into the noise. 

Figure [5] shows the histogram of (2J 7 hl) for the 164971 
surviving candidates up to (2^-"hl) -values of 9. The dis- 
tribution actually extends up to (2J 7 hl) ~ 542.8 but we 
show only the low (2J-HL)-values distribution in order to 
explain our next choice of threshold. There are 166 can- 
didates with (2J 7 hl) > 9 and they are all clustered at 
the two frequencies of ~108.9 Hz and ^575.2 Hz, cor- 
responding to two simulated signals injected in the data 
stream, as discussed in Sec. VI The region (2J"hl) < 6.5 
contains well over 99% of the candidates and, as seen in 
Figure [5] below 6.5 the density of candidates increases 
very sharply. We will take 6.5 as a threshold for the next 
step in our follow-up procedure. 

The number of candidates expected to survive the 6.5 
cut on (2.7-hl) hi fixed 0.5 Hz bands increases with fre- 
quency because the number of sky locations searched 
scales with the square of the searched frequency. In or- 



15 



10 



10 



10 



<2F HL > 10 



<2F H > • <2F L > • <2F HL > 



10 



10 



: 4 




10 



FIG. 4. Values of 2J r averaged over 121 data segments 
for the single-detector case, (2J-h) (green dots), (2J-l) (blue 
dots) and the multi-detector case, (2J-hl) (red dots), against 
those for the combined multi-detector statistic. The top (bot- 
tom) plot shows such values for 164 971 (7 067) surviving (ve- 
toed) S5R5 candidates such that (2J-h) and (or) (2J-l) is less 
(greater) than (2Jhl). 



der to compute the false alarm probability corresponding 
to this threshold in different frequency bands, note that 
in the absence of a signal, the value of 2J 7 hl in the i th 
segment, 2.7^, follows a x 2 distribution with 4 degrees 



HL 

of freedom. Furthermore, since 



121 



(2J- HL ) x 121 = 2^ 



(0 

HL ' 



(41) 



it is clear that (2J 7 hl) x 121 is a x 2 random variable with 
(4 x 121) degrees of freedom. The false alarm probability 
corresponding to a threshold at (6.5 x 121) for such a ran- 
dom variable is ~10 -16 . This corresponds to expected 
false alarm rates of about 0.1%, 0.6%, 2.6% and 10% 
for searches in Gaussian noise over over 0.5 Hz bands at 
100 Hz, 250 Hz, 500 Hz and 1 000 Hz, respectively, con- 
sidering the number of independent trials given by the 
number of searched templates in the respective bands. 
Disregarding the non-Gaussian line features evident in 
Fig. [6j the ratio of the number of candidates observed 
above the 6.5 threshold at lower frequencies (say below 
800 Hz) to that at higher frequencies (say above 800 Hz) 
is not inconsistent|f]with the ratios of the false alarm rates 



4 Due to the low number statistics it is hard to make a sharper 
statement. 



computed above. We note that the false alarm probabil- 
ity given above overestimates the number of expected 
candidates above the 6.5 threshold because that thresh- 
old is not the only cut applied to the data. The previous 
cuts, discussed in the preceding sections, lower the actual 
false alarm probability of the surviving candidates. 
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FIG. 5. Histogram of average multi-detector 2 J-"- values for 
164 971 S5R5 surviving candidates. The red dotted line draws 
the boundaries of the bulk of candidates due to instrumental 
noise, and corresponds to the threshold (2_Fhl) = 6.5. 



There are 184 remaining S5R5 candidates, for which 
(2.Fhl) > 6.5 and they are shown in Fig. [7] as a function 
of the frequency. They are clustered at twelve frequen- 
cies and only the most significant candidate from every 
cluster has been followed up. 



E. Following up candidates with S5R3 data 

For the next step in our post-processing pipeline, recall 
that our underlying signal model of Eq. ([I]) assumes that 
the signal is long-lived; thus its amplitude is constant 
in time and its intrinsic frequency evolves smoothly ac- 
cording to Eq. ([5]) with a constant spindown. This is 
an idealized model: although pulsars are the most sta- 
ble clocks in the Universe, neutron stars are nonetheless 
known to glitch, to be perturbed by external agents, and 
in some cases to be affected by significant timing noise. 
Furthermore, for sufficiently long observation times, the 
spin-down evolution model that we use here, including 
only the first spindown order, may not be adequate to 
describe the actual signal model (see also the discussion 
in Section III D ) . However, since the data set used in the 



S5R3 run ends just about a week before the S5R5 data 
set, it is reasonable to assume that any putative signals 
should be present in both data sets. Moreover, the av- 
erage noise floor level in the detectors turns out to be 
approximately stable between S5R3 and S5R5. Thus, we 
might expect that a detectable signal should be visible in 
both searches. 
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FIG. 6. Average multi-detector 2T- values of the S5R5 candi- 
dates surviving the ./-"-statistic consistency veto as a function 
of the frequency. The horizontal line represents the threshold 
value of (2.7-hl) = 6.5. The bottom plot shows the top plot 
in the region close to the threshold. 



The next follow-up step for each candidate then con- 
sists of a hierarchical search carried out on the same WU 
as done for S5R5 (i.e. over the same parameter space), 
but using the S5R3 data set. The closes^ candidate to 
the original one from such a search was identified and the 
value of its detection statistic was compared with what 
one would expect if the S5R5 candidate were due to a sig- 
nal. In particular, the expected number counts in S5R3 
and S5R5 should be related according to 



£rjy£"W] = x 84/121 



(42) 



where 84 is the number of data segments used in S5R3, as 
we shall see in Appendix [B| Possible reasons for this not 
to be a good approximation would be if the detector noise 
floor were to vary significantly between S5R3 and S5R5, 
or the fact that the relative geometry between the detec- 
tor and source varies in time. We have already remarked 
that the noise floor is, on the average, stable between 
S5R3 and S5R5. Furthermore, albeit in each segment 



5 The distance used to judge closeness between candidates is a 
Euclidean distance expressed in bins in the four dimensions 
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FIG. 7. Histogram of the S5R5 candidates that have been 
further investigated through deep follow-up study with the 
S5R3 data set. 



the expected 2J r -values for a given signal might be dif- 
ferent, if we average this expected value over many non- 
overlapping segments we expect this to converge within a 
few tens of segments; recall here that each segment spans 
a duration of 25 hours, while the antenna pattern func- 
tion of the detectors has a periodicity of 24 hours. Thus, 



it is reasonable to assume that Eq. ( 42 1 is valid. 

Candidates for which the measured value 7V^ 5R3 was 
more than 3 a less significant than the expected £'[A r ( ? 5R3 ] 
were discarded as not being consistent with a CW sig- 
nal, where a was computed using Eq. ( 22 ) . As shown in 
Table [TTJ two candidates were discarded by this follow- 
up test, at ^80.9 Hz and ~108.9 Hz. However, the 
second of these, as well as the candidates at ~52.8 Hz 
and ^575.2 Hz, represent three simulated signals injected 



(/,/,«, S). 



only part of the time during S5, as discussed in Sec. |VI| 
As we can see from Fig. [7J the bulk of the candidates 
arise from the strong hardware-injected pulsar 2 and 3 
signals (see Sec. VI I. In particular, 84 and 88 candidates 
are clustered near the frequencies of these two injected 
signals respectively. 



F. Fully-coherent follow-up 

Excluding the hardware-injected simulated signals, the 
post-processing up to this point has left us with 8 sur- 
viving candidates. These have been significant enough to 
pass our thresholds and have not been clearly identified 
as instrumental artifacts, or eliminated by inconsistency 
between the H and L detectors, or by the follow-up with 
the S5R3 data set. We therefore need to consider other 
more sensitive methods. If these candidates are real sig- 
nals, then their SNRs and significance should increase 
if the parameter space grids are made finer, or as the 
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96.581099597 
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72 
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37 


108.85716501 


13.9 


101 


70 


55 


144.74321811 


7.9 


71 


49 


42 


434.09886421 


7.7 


70 


49 


46 


575.16357663 


11.6 


89 


62 


53 


677.47882796 


8.0 


71 


49 


46 


932.36948703 


8.5 


74 


51 


46 


984.44286823 


8.2 


73 


51 


47 


1030.1650892 


9.1 


77 


53 


53 


1141.9926498 


9.4 


78 


54 


48 



TABLE II. S5R5 candidates followed-up using the S5R3 data 
set. The different columns represent: (I) the candidate fre- 
quency (shown also in Fig.[7|, (II) its significance (CR), (III) 
the S5R5 number count (iVf 5R5 ), (V) the expected and (VI) 
observed number count values after following the candidate 
up (£[iV c S5R3 ] and JV C S5R3 , respectively). 



coherent integration time becomes larger. 

We use a three-step procedure consisting of a grid- 
based semi-coherent Hough search, followed by a semi- 
coherent and a final fully-coherent ^-statistic search, us- 
ing the Mesh Adaptive Direct Search (MADS) algorithm 
for constrained optimization. The reference implementa- 
tion of the MADS algorithm is publicly available through 
the NOMAD library @H1 HZ]- Hence, in the following, 
we refer to such searches simply as NOMAD searches. 
Contrary to the traditional grid-based methods, a mesh 
adaptive search constructs the trial points as the search 
evolves aiming to find the maximum of the statistic. 

The three steps of the follow-up procedure are the fol- 
lowing: 



1 



Re- run the Hough search around a given candidate, 
but with a finer grid to reduce the mismatch with 
a putative signal. The search region includes 5 fre- 
quency bins on either side around the candidate 
and 16 neighboring coarse sky-grid points. The fine 
Hough sky-grid is refined by a factor of 2 in each 
direction by using p = 1 (see Eq. 13 ) instead of 0.5 
as in the original search. Furthermore, we refine 



the coarse / grid spacing of Eq. (|10| by a factor 
N scg = 121 



2. The loudest candidate from the first step is used as 
a starting point for the semi-coherent ^-statistic 
NOMAD optimization. The detection statistic in 
this step is the sum of the ^-statistic values from 
each segment. This search has been performed in 
a fixed parameter space box around the starting 
point. The dimensions of the box are Af = 10~ 4 
Hz, Aa = 0.10 rad, AS = 0.24 rad and Af = lO" 10 
Hz s _1 . The loudest candidate found in this semi- 
coherent ^-statistic NOMAD search is passed on 
to the next step. 

3. In the third step, the loudest candidate from the 



previous step is used as a starting point for the 
fully-coherent ^-statistic NOMAD search. This 
search spans the entire duration of the S5R5 data 
set and has been carried out in a parameter space 
box denned by using the diagonal elements of the 
inverse Fisher matrix in each dimension around the 
starting point. These elements are described by 
Eq. (16) in |48j and have been computed from the 
inverse of the semi-coherent parameter space metric 
computed at the candidate point, re-scaled by the 
measured SNR at the same point; for more details 
we refer the reader to [4"8] , 

In both the semi-coherent and fully-coherent NOMAD 
searches, we ran multiple instances of the algorithm iter- 
ating over the mesh coarsening exponent using both de- 
terministic [35] and stochastic [35] methods for the choice 
of search directions. Based on Monte-Carlo studies, the 
false-dismissal probability of the follow-up procedure is 
found to be less than 10%. 

As said earlier, in the presence of a real signal we ex- 
pect the significance of a candidate to increase as the tem- 
plate grid becomes finer because there will be a template 
with a smaller mismatch with respect to the real signal. 
At that template the signature of the signal should be 
more evident and all the consistency tests should con- 
tinue to hold. If the candidate signal detected on the 
finer grid does not pass a consistency test, this indicates 
that it is not behaving as we would expect from the sig- 
nals that we are targeting. The candidates at —96.6 Hz, 
-144.7 Hz, -932.4 Hz, -1030.2 Hz and -1142 Hz fail 
a multi-detector versus single-detector ^-statistic con- 
sistency test (see Sec. IV C I after performing the semi- 
coherent NOMAD search; therefore, they cannot be con- 
sidered defensible CW signals. Moreover, line artifacts 
appear in the average power spectrum of S5 H data at 
-932.4 Hz, -1030.2 Hz and -1142 Hz. 

The remaining candidates, namely at —434.1 Hz, 
—677.5 Hz and —984.4 Hz, survive the ^"-statistic consis- 
tency test on the finer grid and are followed-up with the 
fully-coherent ^"-statistic NOMAD search. However, for 
each of them, the maximum value of the detection statis- 
tic over the parameter space searched is much lower than 
would be expected based on the original candidate pa- 
rameters, and in fact is consistent with the expectation 
for Gaussian noise. Hence, also these three candidates 
do not survive a more sensitive inspection and cannot 
be considered viable detection candidates. Thus, we see 
that all the candidates listed in Table |n] are inconsistent 
with the properties of a true CW signal. 



V. UPPER LIMIT ESTIMATION 

The analysis of the Einstein@Home searches pre- 
sented here has not identified any convincing CW sig- 
nal. Hence, we proceed to set upper limits on the max- 
imum intrinsic gravitational wave strain ho that is con- 
sistent with our observations for a population of CW 
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signals described by Eq. from random positions 

in the sky, in the gravitational wave frequency range 
[50.5, 1 190] Hz, and with spindown values in the range of 
~ [—20, 1.1] x 10~ 10 Hz s _1 . The nuisance parameters 
cos I, 4>q, and ip are assumed to be uniformly distributed. 
As commonly done in all-sky, all-frequency searches, the 
upper limits are given in different frequency sub-bands 
and here we have chosen these to be 0.5 Hz wide. Each 
upper limit is based on the most significant event from 
the S5R5 search in its 0.5 Hz band. 



A. Monte-Carlo upper-limit estimates 

Our procedure for setting upper limits uses Monte- 
Carlo signal injection studies using the same search and 
post-processing pipeline (except for the S5R3 and fully- 
coherent follow-ups) that we have described above. In 
every 0.5 Hz band, our goal is to find the value of ho (de- 
noted /Iq 0% ) such that 90% of the signal injections at this 
amplitude would be recovered by our search and are more 
significant than the most significant candidate from our 
actual search in that band. We can thus exclude, with 
90% confidence, the existence of sources (from our spe- 
cific population) that have an amplitude ho > h\ w% . 

In order to estimate /iq 0% , for each injection at a 
randomly chosen parameter space point, a hierarchical 
search is performed over a small parameter space region, 
which consists of: 

• a 0.8 mHz frequency band centered at the S5R5 
frequency grid-point closest to the randomly chosen 
source frequency; 

• 4 spindown values around the S5R5 frequency 
derivative grid-point closest to the randomly cho- 
sen spindown; 

• a sky-patch consisting of 10 S5R5 coarse sky-grid- 
points closest to the randomly chosen sky location. 

At the end of each hierarchical search, the most sig- 
nificant candidate is selected and post-processed as de- 
scribed in the previous section. The vetoes for ex- 
cluding known instrumental lines are not required be- 
cause we have already excluded them from the upper- 
limit analysis. We first check if this candidate is sig- 
nificant enough to be part of the 100 Einstcin@Home 
loudest candidates originally selected in its correspond- 
ing 0.5 Hz band. Then we perform the ^-statistic consis- 
tency check and finally we compare the computed aver- 
age multi-interferometer 2 J 7 - value ((2-7 r §f md )) with the 
maximum (2J 7 hl) value we have in the corresponding 
0.5 Hz band. If (2J r §£ nd ) is greater than the maximum 
(2J"hl), then the simulated source is considered to be 
recovered and more significant than the most significant 
candidate of the search. The confidence level is defined 
as C — n ICC /n to t, where n rcc is the number of recovered 
candidates, and n tot = 100 is the total number of injec- 
tions performed. 



After some preliminary tuning to determine a range 
of h values close to the 90% confidence level, we use 
an iterative procedure to determine the confidence as a 
function of the injected population until we hit a con- 
fidence value close to 90%, within the expected 1 a fluc- 
tuations. Since we use 100 injections, from a binomial 
statistic we estimate the 1 a fluctuation to be 3% and 
hence we associate the ft,Q 0% value to any measured con- 
fidence in the range 87% — 93%. The 3% uncertainty in 
confidence translates to an uncertainty in /iq ^ smaller 
than 5%, as can be seen from Fig. [8| which shows a typ- 
ical confidence versus injected behavior. Each point 
in Fig. [8] was derived with 1 000 injections and hence is 
affected by fluctuations smaller than 1%. 
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FIG. 8. Confidence versus the injected ho values for sets of 
1000 injections in the band [216,216.5] Hz. This plot illus- 
trates how the uncertainty on the confidence level affects the 
uncertainty on the value of hg 0/ °. 

The lower (red) curve in Fig. [9] shows the resulting 
upper limits as a function of the frequency. The upper 
(blue) curve shows the upper limit values from the previ- 
ous Einstein@Home search in early S5 data [24 . The 
current upper limit values are about a factor 3 more 
constraining than the previous Einstein@Home ones. In 
particular, the most constraining upper limit falls in the 
0.5 Hz-wide band at 152.5 Hz, where we can exclude the 
presence of signals with h greater than 7.6 x 10~ 25 . The 
three stars shown in Fig. [9] correspond to the simulated 
pulsars 2, 3 and 5, i.e. the hardware injections recovered 
in the S5R5 search (discussed in Sec. |VI[). 

The numerical data for the plot in Fig. [9] can be ob- 
tained separately ^\ A conservative estimate of the over- 
all uncertainty on the /Iq 0% values shown in Fig.Nis 15%, 



See Supplemental Material at [Replace PRD link instead than 
ID LIGO-P1200026] for numerical values of upper limits. 
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FIG. 9. Upper limits for the S5R5 EinsteinOHome search (red dots) as well as the previous EinsteinQHome search, called 
S5R1, which used early S5 data (blue dots) |24j . The three stars correspond to hardware-injected simulated pulsars which 
were recovered in the S5R5 search. The curves represent the source strain amplitude ho at which 90% of simulated signals 
would be detected. The vertical bars represent 156 half-Hz frequency bands contaminated by instrumental disturbances for 
which no upper limits are provided. The upper limits for the 0.5 Hz-wide bands starting at 69.5 Hz, 139.5 Hz, 335.5 Hz and 
648 Hz are fairly high due to significant partial contamination in these bands by lines listed in Tables |VI| and |VII| Note that 
the broadness of the red curve is due to the 5% steps used to vary the injected population ho values in the Monte-Carlo signal 
software-injections until a confidence value close to 90% is reached. In addition, less than 1 /4 of the spectral range shown was 
excluded in many narrow bands because of kn own instrumental artifacts, as described in Sec. |IVB| The cyan curve shows the 
predicted ho° % upper limits according to Eq. (431. 



having added to the 1 a statistical upper limit estimation 
procedure uncertainty the 10% amplitude calibration 
uncertainties for the data used in this Einstein@Home 
run [50]. 

As we have excluded from the search those frequency 
bands hosting spectral artifacts (Tables VI and VII) 
and the cleaned noise bands (Table 0, 

we therefore also 

exclude these frequency intervals from the upper limit 
statements. Vertical bars in Fig. [9] represent 156 half-Hz 
frequency bands for which no upper limits are provided 
because the entire half-Hz band has been excluded. 

As shown in Fig. [9j the upper limits on h provided 
in the 0.5 Hz-wide bands starting at 69.5 Hz, 139.5 Hz, 
335.5 Hz and 648 Hz are fairly high, roughly equal to 3 x 



10 



-2.S 



This is due to significant partial contamination in these 
bands by lines listed in Tables |VI| and |VII| the upper 
limit is given for the remaining, clean part of the band, 
but loud candidates from the disturbed part make up 
the loudest 100 candidates selected in the processing, so a 
simulated signal must be especially loud to surpass those. 
Note that, for the same reason, if we had set upper limits 
in the 156 half-Hz bands shown in Fig. |9]we would have 
obtained similarly high upper limits on h - 



B. Analytic sensitivity estimates 



The hn ^ upper limits can be independently predicted 



:10 24 , 1.7x10 23 and 2x10 23 , respectively. using the method in [5T], adapted to the Hough-on- J 7 - 
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statistic search method (see [55] for details). The upper 
limit procedure described above is modelled by a sim- 
ple threshold on the number count, where the thresh- 
olds are the largest number counts observed in each 
0.5 Hz upper limit band. The probability that, in the 
neighbourhood of an injected signal, the number count 
n c will exceed a threshold n Cjt h is denoted by P[n c > 
n c ,th\p(ho, cos l, 4>o, ip, Tn)]; this probability can be calcu- 
lated analytically from the known distribution of n c . The 
recovered SNR, p, is a function of the nuisance parame- 
ters and of the mismatch m between the injected signal 
and the nearest template. In addition, note that in the 
presence of a signal, 2J 7 follows a non-central \ 2 distri- 
bution with 4 degrees of freedom; p 2 is the non-centrality 
parameter of this distribution. In the presence of a sig- 
nal, averaging P over the parameters of p (except ho) 
gives (P(n c > n Cl tb\p(ha))) ■> which equals the confidence 
of recovering a population of injections with amplitude 
ho. In each 0.5 Hz band, we determine the value of ho 
such that (P{n c > n c ,th\p(ho))) = 90%; this value is then 
the predicted value of h w% , given by 

h 9 0% = (43) 

V J data 

as a function of the detector noise Sh and the total data 
volume Tdata = N scg T scg . The factor H varies from ^141 
to ~ 150 over the range of search frequencies, and is plot- 
ted in [52 . It is given by H = 2.5p 90% y^N scg , where 
p 90% is the mean injected SNR per segment of a popula- 
tion of signals as described above, and is itself a function 
of the false alarm and false dismissal probabilities, and 
N seg jHIl [52] . The variation of H as a function of fre- 
quency arises from the variation of p 90% function 
of the false alarm probability in each upper limit band, 
which are calculated from the largest number counts n c 
plotted in Fig. [3] 

The predicted values are shown by the cyan curve in 
Fig. [9] The root-mean-square error between the Monte- 
Carlo estimated and predicted /ig 0% values (~ 7% over 
all frequencies) is comparable with the uncertainties due 
to calibration (10%) and the finite stepping of the Monte 
Carlo procedure (5%). This demonstrates that the sen- 
sitivity of the Hough-on- ^-statistic search method, as a 
function of the search parameters, is well understood. 
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FIG. 10. Panel (a) and (b) represent the distance range (in 
kpc) and the maximum ellipticity, respectively, as a function 
of the frequency. Both the panels are valid for neutron stars 
spinning down solely due to gravitational radiation and as- 
suming a spindown value of ~ — 2 nHz s _1 . In these plots, 
the 156 half-Hz frequency bands for which no upper limits are 
provided have not been considered. 



derivative /, and emitting at the spindown limit is 



C. Astrophysical reach 

Figure [10] shows the maximum reach of our search. 
The top panel shows the maximum distance at which we 
could have detected a source emitting a CW signal with 
strain amplitude /jq 0% . The source is assumed to be spin- 
ning down at the maximum rate considered in the search, 
~ —2 nHz s _1 , and emitting at the spindown limit, i.e. 
with all of the lost rotational energy going into grav- 
itational waves. The intrinsic gravitational wave strain 
from a source at a distance d, with frequency /, frequency 



hf = 4.54 x 10 



-24 



'1 kpc\ / 250 yr \ 1/2 

~) W/(W 



(44) 



where the canonical value of 10 38 kg m 2 is assumed for 
I zz in Eq. Q. The bottom panel of the figure does not 
depend on any result from the search. It shows the spin- 
down ellipticity values as a function of the frequency for 
sources emitting in gravitational waves all the energy lost 
while spinning down at a rate of ~ — 2 nHz s _1 . This 
is obtained by setting / = —2 nHz s _1 in the following 
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equation: 



e sd = 7.63 x 10" 



Pulsar 



Pulsar 1 



10 



-10 



Hz s- 



/ 100 Hz^ 5/ ~ 
/ 



(45) 

Around the frequency of greatest sensitivity, 152.5 Hz, 
we are sensitive to objects as far as 3.8 kpc and with an 
ellipticity e ~ 10 -4 . Normal neutron stars are expected 
to have e less than a few times 10 -6 based on theoretical 
predictions [31]. A plausible value of e ~ 3.5 x 10~ 6 
could be detectable by a search like this if the object 
were emitting at 625 Hz and at a distance no further 
than 500 pc. 



VI. STUDY OF HARDWARE-INJECTED 
SIGNALS 

As part of the testing and validation of search pipelines 
and analysis codes, simulated signals are added into the 
interferometer length control system to produce mirror 
motions similar to what would be generated if a gravi- 
tational wave signal were present. Table |III| shows the 
parameters of the set of simulated CW signals injected 
into the LIGO detectors; we shall often refer to these in- 
jections also as "fake pulsars" . These injections were ac- 
tive from the GPS epoch 829 412 600 s until 875 301 345 s. 
Of these ten hardware-injected CW signals, eight had fre- 
quencies covered by the S5R5 search frequency band: the 
fake pulsars 4 and 7 have frequencies outside this band 
and thus have not been taken into account during this 
analysis. The fake pulsars 6 and 8 have spindown values 
outside the S5R5 search frequency derivative range. 

As a minor complication, the hardware injections were 
not active all the time. In the S5R5 data set, their duty 
cycle was ~ 63% and ~ 60% in L and H, respectively. 
The hardware injections were active in 76 of the 121 
S5R5 data segments, and they were completely absent 
in the remaining 45 segments. Their expected value of 
the Hough number count is thus: 



E\n 



HIi 



3 



i 



(46) 



where the superscript "HI" refers to Hardware Injection, 
j runs over the number of data segments where the hard- 
ware injections were active and I runs over the remaining 
segments, and the Wj are the Hough weights given by 



Eq. (18); r]j and p are the probabilities that the esti- 
matedvalue of 2J 7 (for a given data segment) crosses the 
threshold 2J r t h = 5.2 in the presence and absence of a sig- 
nal, respectively (expressions for r]j and p can be found 
in [EJ [T§]). 

Figure [IT] shows, for each pulsar hardware injection, 
the histograms of the expected 2J r -values, computed for 
all the segments where a particular hardware injection 
was active. From this figure one can infer that the fake 
pulsars expected to be the loudest in terms of 2J r -values 
are pulsars 2, 3, 5 and 8. 
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FIG. 11. Histograms of the 1J-- values estimated, for every 
data segment, for all the S5 hardware injections covered by 
the investigated frequency range. The largest 2J r -values are 
from the fake pulsars 2, 3, 5 and 8. In all the plots, the total 
number of data segments where the hardware injections were 
active is 76. 



The pulsar hardware injections go through the normal 
search and post-processing pipelines in the usual way as 
described in the previous sections. As discussed in the 
following, the fake pulsars recovered in the S5R5 search 
are pulsars 2, 3 and 5. On the other hand, fake pulsars 
6 and 8 are missed because their spindown values are 
outside our search range. As expected, fake pulsars 0, 
1 and 9 are not recovered because their amplitudes are 
too weak and they do not pass the (2J 7 hl) > 6.5 cut, as 
shown in Table |IV| In this table, the observed number 
counts of 69 to 72 are consistent with noise fluctuations 
in those half-Hz bands (see Fig. [3]) superseding the weak 
injected signals. Note that the expected 2J r -values shown 
in Fig. [TT] were computed assuming a search using exactly 
the correct signal parameters provided in Table [HT[ while 
the observed (2J 7 hl) -values in Table IV were obtained 



from our search which uses a grid of templates, so signif- 
icant mismatch is to be expected. Table [TV] compares the 
expected and observed values of the number counts asso- 
ciated with the S5 hardware injections and the surviving 
S5R5 candidate events closest to them. As in Sec. |IV E 



the measure of distance used here is a Euclidean distance, 
expressed in bins, in the four dimensions (/, /, a, S). 

Table [V] lists the parameters of pulsars 2, 3 and 5 
and the parameters of the corresponding recovered can- 
didates. We successfully find candidates near the cor- 
rect signal parameters. The mismatch in spindown might 
seem large, but in fact the injections were found at the 
nearest spindown template. The number count values 
show consistency within the 3er range. 
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Name 


fp (Hz) 


/ (Hz s" 


l ) 


a (rad) 


5 (rad) 


ip (rad) 


h (rad) 


cos l (rad) 


h 







Fake pulsar 


265.5771052 


-4.15 x 10" 


1 2 


1.248817 


-0.981180 


0.770087 


2.66 


0.794905 


2.47 x 


10" 


2 b 


Fake pulsar 1 


849.0832962 


-3.00 x 10" 


10 


0.652646 


-0.514042 


0.356036 


1.28 


0.463822 


1.06 x 


10" 


24 


Fake pulsar 2 


575.163573 


-1.37 x 10" 


13 


3.756929 


0.060109 


-0.221788 


4.03 


-0.928576 


4.02 x 


10" 


24 


Fake pulsar 3 


108.8571594 


-1.46 x 10" 


17 


3.113189 


-0.583579 


0.444280 


5.53 


-0.080666 


1.63 x 


10" 


-23 


Fake pulsar 4 


1403.163331 


-2.54 x 10 


—8 


4.886707 


-0.217584 


-0.647939 


4.83 


0.277321 


4.56 x 


10" 


-23 


Fake pulsar 5 


OZ.oUooZ4oO 


- — 4.UO X 1U 


18 


O.Zolool 


— i.40ozoy 


— u.oDoyoo 


Z.Zo 




4.O0 X 




24 


Fake pulsar 6 


148.7190257 


-6.73 x 10 


-9 


6.261385 


-1.141840 


0.470985 


0.97 


-0.153733 


6.92 x 


10" 


-25 


Fake pulsar 7 


1220.979581 


-1.12 x 10 


-9 


3.899513 


-0.356931 


0.512323 


5.25 


0.756814 


2.20 x 


10" 


24 


Fake pulsar 8 


194.3083185 


-8.65 x 10 


-9 


6.132905 


-0.583263 


0.170471 


5.89 


0.073903 


1.59 x 


10" 


-23 


Fake pulsar 9 


763.8473165 


-1.45 x 10" 


17 


3.471208 


1.321033 


-0.008560 


1.01 


-0.619187 


8.13 x 


10" 


-25 



TABLE III. Simulated ("fake") pulsar hardware injections during the S5 LIGO run, created with the JPL DE405 Sun and 
Earth ephemeris files. The pulsar parameters are defined at the GPS reference time of 751680 013 s in the SSB frame. These 
are the frequency fp, the spindown /, the sky position (a, 5), the polarization angle tp, the initial phase (j>o, the inclination 
parameter cost and the dimensionless strain amplitude ho- These parameters correspond to the only set of hardware injections, 
injected into the S5 LIGO data, that fall within the GPS times of the S5R5 data. 



TABLE IV. Comparison between the expected and observed 
number counts (E[n^ 1 ] and N^ 5RS respectively) associated 
with the hardware injections and the recovered S5R5 candi- 
dates closest to these injections. The (2J-HL)-values for each 
of the candidates are also listed. The fake signals labeled 
pulsar 4 and 7 were not taken into account in this analysis 
since they have frequencies outside the S5R5 search frequency 
range. The "expected" values marked by asterisks are not ac- 
tually expected to be obtained because the spindown rates for 
those signals lie outside the range of this search. 



Name 


E[n™] 


iVf R5 


(ZFhl) 


Fake pulsar 


53 


69 


5.6 


Fake pulsar 1 


49 


71 


5.7 


Fake pulsar 2 


88 


80 


22.3 


Fake pulsar 3 


89 


96 


197.1 


Fake pulsar 5 


85 


70 


6.6 


Fake pulsar 6 


76* 


71 


5.3 


Fake pulsar 8 


87* 


72 


5.5 


Fake pulsar 9 


53 


72 


5.8 



TABLE V. Study of hardware injections in the S5R5 search. 
Values of (/, a, S, /) for the fake pulsar 2, 3, 5 and for the clos- 
est recovered candidates (denoted as Cand 2, 3, 5) are listed. 
The n c -value in the last column represents the expectation 
value -Bfn^ 1 ] for the fake pulsars and the observed number 
count iV,? 5 for the corresponding recovered candidate. 



Name 




/ (Hz) 


a (rad) 


<5 (rad) 


/ (Hz s" 1 


) 




Fake pulsar 
Cand 2 


2 


575.163573 
575.163556 


3.756929 
3.757514 


0.060109 
0.065354 


-1.37 x 10~ 
-1.64 x 10" 


-13 
-11 


00 00 
O 00 


Fake pulsar 
Cand 3 


3 


108.857159 
108.857158 


3.113189 
3.09806 


-0.583579 
-0.5839483 


-1.46 x 10" 
-1.64 x 10" 


-17 
-11 


89 
96 


Fake pulsar 
Cand 5 


5 


52.8083243 
52.8082977 


5.281831 
5.58845 


-1.463269 
-1.470972 


-4.03 x 10" 
-1.64 x 10" 


-18 
-11 


85 
70 



VII. CONCLUSIONS 

No evidence for continuous gravitational waves has 
been observed in the search presented here. Upper lim- 
its on the intrinsic gravitational wave strain have been 
derived using standard population-based methods, and 
are shown in Fig. [9] These results are about a factor of 
3 more constraining than those from the previous Ein- 
stein@Home search in early S5 data [24] . This improve- 
ment comes from using more data (a year versus two 
months), from using a multi-detector coherent statistic 
(versus a single-detector statistic), from a lower effec- 
tive threshold in the coherent detection stage, and from 
a more sensitive incoherent method to combine the in- 
formation from the coherently analyzed segments. The 
largest effect comes from lowering the effective thresh- 
old. Indeed much of the improvement in sensitivity is 
attributable to improved data analysis methods (as op- 
posed to improved detector sensitivity). If we had used 



the much higher threshold of 25 on 2J 7 , as in [23], our 
sensitivity would have been a factor of ~ 2.5 worse than 
our final upper limits, thereby undoing almost all of the 
factor of 3 improvement mentioned above. 

We have not included second time-derivatives of the 
frequency in our search. This could be astrophysically 
significant in some regions of parameter space, as dis- 



cussed in Sec. Ill D It is important to keep this caveat 
in mind while interpreting our results. 

This is the most sensitive wide-frequency-range, all- 
sky search for CW signals performed to date. The upper 
limit values are comparable to those obtained recently 
using the PowerFlux method [22J [2S] on the entire S5 
data set (S5R3+S5R5). Ref. [26] searched for CW sig- 
nals over the whole sky, in a smaller frequency band (up 
to 800 Hz versus 1 190 Hz here), but a broader spindown 
range up to —6 nHz/s. Strain upper limits were set at the 
95% confidence in 0.25-Hz wide sub-bands. In particular, 
near 152 Hz, the PowerFlux strict, all-sky upper limit on 
worst-case linearly (best-case circularly) polarized strain 



23 



amplitude ho is ~ 1 X 10 24 (3.5 x 10 25 ). As a com- 
parison, at the same frequency, this search constrains the 
strain to h < 7.6x 10~ 25 (as shown in Fig.^, 9.2x 10~ 25 
and 3.2 x 10 -25 for the case of average, linear and circular 
polarization, respectively, with a 90% confidence level in 
a 0.5-Hz wide band. 

The most constraining upper limit obtained by this 
search is /ig 0% ~ 7.6 x 10~ 25 at 152 Hz; the corresponding 
maximum reach is roughly 4 kpc, assuming e ~ 10~ 4 . 

It has long been expected that searching a large pa- 
rameter space for CW signals will require hierarchical 
semi-coherent searches. This analysis is a milestone to- 
wards that goal, and we expect that future analyses will 
build on the tools developed here. 

There are a number of areas where further improve- 
ments are possible. In the latest round of analysis 
(an Einstein@Home processing run that began in March 
2012), some of the post-processing techniques developed 
for this analysis have been "moved upstream" to the 
hosts. One example is the generalized /"-statistic con- 
sistency test [53] . This continues the pattern of moving 
analyses formerly carried out in the post-processing stage 
(for example, the incoherent combination step) onto host 
machines. Another step forward is in the semi-coherent 
algorithm that combines the coherent analyses from the 
different segments. The Hough algorithm described here 
turned out to be rather cumbersome, and does not com- 
bine the coarse and fine grids in an optimal way. The 
latest round of Einstein@Home processing makes use of 
a simpler optimal semi-coherent method, which allows 
longer coherent time baselines to be used. This method, 
based on a detailed analysis of correlations in parameter 
space [51], is described inJSSJ. Looking farther forward, 
we expect to use higher-order spindown parameters both 
to search for a broader class of signals as well as to be able 
to employ longer coherent time baselines in the analysis. 

The Advanced LIGO and Advanced Virgo detectors 
are currently under construction, and should begin op- 
erations around 2016. In comparison with the current 
generation, these instruments will provide an order-of- 
magnitude improvement in strain sensitivity, increasing 
the volume of space observed by a factor of a thousand. 

These and other improvements in data analysis meth- 
ods and instrumentation make us optimistic that we will 
eventually be able to make direct detections of CW sig- 
nals. Such detection will provide new insights into the in- 
ternal structure, formation history and population statis- 
tics of neutron stars. 
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Appendix A: Problems in calculating the weights 

In this section we describe two software errors that 
affected the main hierarchical search code. While these 
errors did not invalidate the search results and were dealt 
with adequately, we document them here for complete- 
ness. 

The first issue is connected with the choice of weights 



used to construct the number count defined in Eqs. (17) 
and (18). The weights for each segment are computed 



following Eq. ( 19 ) which uses the harmonic mean of the 



noise spectra for each SFT. However, the original method 
for computing the weights used the arithmetic mean of 
the contributions instead of the harmonic mean, which 
turned out to have the perverse effect that segments with 
a few noisy SFTs got a disproportionately large weight. 
This led to a much larger value for the variance given in 
Eq. (22) and a correspondingly anomalously small value 
of the significance defined in Eq. (21). All WUs were 



originally run using the arithmetic mean which led to 
anomalously low values of the significance CR due to 
non-stationary noise for a small number of frequency 
bands: [50, 202] Hz, [328.5, 329] Hz, [995.5, 1 010] Hz and 
[1069,1075] Hz. The WUs for these bands were re- 
computed with the weights given by Eq. ( |19| ). All other 
frequencies are unaffected by this problem. 

A second issue, which interacts with the problematic 
calculation of the weights described above, is floating 
point inaccuracy in our implementation of the Hough- 
transform algorithm. A single threshold crossing for the 
^-statistic leads to a +1 in number count for possibly 
a large number of points in parameter space [19] , and it 
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is not necessary to step through parameter space point- 
by-point to calculate the final number count. For the 
vast majority of cases, our implementation of the Hough- 
transform agrees with the brute force approach for calcu- 
lating the number count, but the two can occasionally dif- 
fer. If we were not using weights, these differences would 
have a minor effect on the number count. Occasionally, 
however, these floating point errors coincide with the 
cases when we assign excessively large weights to partic- 
ular segments as discussed above. In these cases, the dis- 
crepancies in the number count can be large and in some 
cases may even exceed the number of segments, which is 
in principle a strict upper bound on the number count. 
Note that our upper limits remain valid because they con- 
sistently use the same search code, and any candidates 
are followed-up by independent codes thereby avoiding 
spurious false alarms. Using the modified weights based 
on Eq. ( 19 ) fixes this problem as well. 



Appendix B: S5R3 Post-Processing 
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FIG. 12. Significance of 230 000 loudest candidates selected 
in 0.5 Hz- wide frequency bands as a function of the frequency. 



The S5R3 run was launched on September 23, 2007 
and ended on September 25, 2008. Like S5R5, it was 
an all-sky search. It used 7 237 S5 LIGO SFTs, col- 
lected between the GPS times of 818 845 553 s (Sat Dec 
17 09:05:40 GMT 2005) and 851 765 191 s (Tue Jan 02 
09:26:17 GMT 2007). The data analyzed consisted of 
3 781 SFTs from H and 3 456 SFTs from L. The num- 
ber of data segments used for the S5R3 run was 84, 
with duration each T scg = 25 hours. The search fre- 
quency range was [50, 1 200] Hz, with a frequency resolu- 
tion Sf ~ 6.7 x 10 -6 Hz and a mismatch value m = 0.3, 
leading to the spindown resolution given by Eq. (16). In 



the S5R3 run, each WU analyzed a constant frequency 
bandwidth B ~ 16 mHz, the full spindown interval, rang- 
ing roughly from —1.6 nHz s _1 to —3.1 x 10~ n Hz s _1 
in steps of 3.8 x 10~ 10 Hz s" 1 , and a certain region of the 



sky, as already discussed in Sec. Ill C The sky-grids and 
output file formats are identical to those used in S5R5. A 
major difference with S5R5 are the weights: no weighting 
scheme was use d in S5R3 and all the weights it), appear- 
ing in Eq. (171 were set to unity. In total, ~ 7 x 10 10 
S5R3 candidates were sent back to the Einstein@Home 
server. The S5R3 post-processing was performed before 
that for S5R5, and it was used as a "test-bed" for the 
latter. It consists of the same series of steps described 



previously in Sec. |IV| The significance values of the 100 
loudest candidates in 0.5 Hz- wide frequency bands are 
plotted in Fig. [12] as a function of the frequency. These 
represent a total of 230 000 candidates. This set is re- 
duced by ~ 27% via the removal of instrumental noise 



artifacts (listed in Tables |VI| and VII ) and of the can- 
didates from search frequency bands close to the fake 
noise (according to Table |l|). For all the surviving 167 779 
candidates, the ^-statistic consistency veto described in 
Sec. IV C removes an additional 3.6%. Figure [T~3] shows 
the values of (2J 7 h), (2Jl) and (2.7-hl) versus (2J 7 hl) 



for 161 785 candidate events that survived (top plot) and 
5 994 that were excluded (bottom plot). 
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FIG. 13. Values of 2T averaged over 84 data segments 
for the single-detector case, (2J-h) (green dots), (2.7%) (blue 
dots), the multi-detector case ((2.7-hl), red dots) against those 
for the combined multi-detector statistic ((2J-hl}). The top 
(bottom) plot shows such values for 161 785 (5 994) surviving 
(vetoed) candidates such that (2J-h) and (or) (2J-"l) is less 
(greater) than (2Jhl}- 



Candidates whose (2.7-hl) -value is greater than 6.5, i.e. 
1 465 out of 161 785, are followed-up by performing a hier- 
archical search using the S5R5 data set. This run, whose 
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details are described in Section IV consists of ~ 46% 
more data than was used for S5R3 and is thus more sen- 
sitive than S5R5. The followed-up candidates are plotted 
in Fig. [14] as a function of the frequency: 87 candidates 
are clustered at ~ 108.857 Hz and represent the contri- 
bution of the hardware injection pulsar 3 [56]; 73 can- 
didates have frequencies peaked around ~ 1081 Hz, 59 
candidates at — 1198.9 Hz and 28 at ~ 1141 Hz. Except 
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FIG. 14. Histogram of 1 465 S5R3 candidates that have been 
further investigated through follow-up study using the S5R5 
data set. 



for the candidates that represent the contribution from 
the hardware injected signal, none of the S5R5 candi- 
dates is consistent with a signal at a level consistent with 
the observed S5R3 excess. We then conclude that no 
gravitational wave signal is observed in the S5R3 data. 



Appendix C: Instrumental noise artifacts 

This appendix contains lists of the main known spec- 
tral lines of instrumental origin in the LIGO detectors 
during the S5 run. They were individually identified with 
a particular source, or were members of identified combs 
of lines found in many channels (e.g. 60 Hz combs), or 
the source was unknown, but they were found in fre- 
quency coincidence between the gravitational wave chan- 
nel and auxiliary channels. In the latter case, to ensure 
that actual gravitational wave signals were not rejected, 
the SNR in the auxiliary channel had to be at least 5 
times larger than in the gravitational wave channel and 
the density of lines in the auxiliary channel had to be 
low enough that accidental coincidence with a line in the 
gravitational wave channel was highly unlikely. The spec- 
tral lines and harmonics detected in H and L are listed 
in Tables |VI| and |VII| respectively. As mentioned ear- 
lier, on the basis of these tables, about 22% and 25% of 
candidates have been excluded from the analysis in the 



S5R5 and S5R3 post-processing, respectively. 

For each candidate with frequency f c and derivative 
f c , the candidate was rejected if a band A/ c ~ f c x 
10~ 4 + |/c| x 10 7 s on cither side of the signal had any 
overlap with an instrumental line band. This was done in 
order to take into account the maximum possible Doppler 
shift due to the Earth's orbital velocity, which is roughly 
10~ 4 in units of the speed of light, and the maximum 
frequency shift due to the spindown of the source over 
the ~ ±10 7 s time span relative to the reference time 
during the EinsteinQHome run. 
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TABLE VI. Known S5 H spectral artifacts. The different 
columns represent: (I) the central frequency of the instru- 
mental line; (II) the number of harmonics including the fun- 
damental; (III) Low-Frequency (LF) bound of the knockout 
band; (IV) High- Frequency (HF) bound of the knockout band; 
(V) the cause of the line (see key below). When there were 
higher harmonics, the third and fourth columns were multi- 
plied by the harmonic number to yield the proper LF and HF 
bounds. 
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Unknown: a line of unknown source that appeared in 
auxiliary channels and met the rejection criteria 
noted in the text. 



TABLE VII. Known S5 L spectral artifacts. The columns are 
the same as in Table IVT1 

A short explanation of the key to the line sources listed 
in the fifth column of Tables IVII and IVHI follows. 

Mains: lines at multiples of the 60 Hz electrical power 
system frequency; the dominant coupling mecha- 
nism at 60 Hz was from magnetic fields generated 
by electric currents coupling to the permanent mag- 
nets mounted on the test masses. 

Electronics: produced by either electronic circuit oscil- 
lations, or by slight data corruption associated with 
repetitive processes in the data acquisition comput- 
ers. The line was identified in power supply voltage 
variation, magnetic fields from electronics or by di- 
rect measurements. 

Calibration: produced for calibration purposes by mov- 
ing a test mass with the actuation system. 

Timing: introduced by the timing verification system. 

Wire: a vibrational resonance of a mirror suspension 
wire. 

Pulsed-heating: produced by cyclically pulsed mains 
heating circuits, coupling to the test mass magnets 
via magnetic fields. 
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